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Pin addition to trajectory calculations, TRAID also provides assistance in the 

following areas: (1) card-data input, (2) printed ¢::tput, with standardized 
formats, (3) vector- and matrix-manipulating routines, (4) problem 
superviston -- e.g., printing title page, enforcing time and page limits, and 
(5) miscellaneous aids -- e.g., plotting, manipulating tabular data. 


If you have a trajectory problem, TRAID can make your programming task easter 
in two ways. First, your program will accompl{sh more per FORTRAN statement. 
The second (hidden) benefit is that your program will probably contain fewer 
errors, because your attention can remain concentrated on the problem as 
distinct from the “spelling.”-© 


4 
To use TRAID, you must write a main program in FORTRAN. In this program you 
must worry about two things: (1) the allocation of storage for arrays, and 
(2) calling the appropriate TRAID routines to perform the desired operations. 
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. 1.0 TRAID -- WHAT AND WHY? 
TRAID is intended to help solve trajectory problems. 
TRAID is a family of subroutines which handle the calculation of 
‘ powered and guided trajectories and of Keplerian orbits. There are sub- 
routines which integrate vehicles forward or backward in time, with thrust- 
mass~-lift-drag-response input by the user, with either preassigned or 
computed guidance commands. Orbits may be elliptical (multiple revolutions 
: are allowed) or hyperbolic; transfer between integrated flight and orbital 
representation is accomplished automatically. 
In addition to trajectory calculations, TRAID also provides assist- 
ance in the following areas: 
card~data input 
printed output, with standardized formats 
vector- and matrix~manipulating routines 
problem supervision -- e.g., printing title page, 
enforcing time and page limits 
miscellaneous aids -- e.g., plotting, manipulating 
tabular data. 
If you have a trajectory problem, TRAID can make your programming 
task easier in two ways. First, your program will accomplish more per 
FORTRAN statement. The second (hidden) benefit is that your program will 
probably contain fewer errors, because your attention can remain concen- 
trated on the problem as distinct from the "spelling." 
¢. 
To use TRAID, you must write a main program in FORTRAN. In this 
' program you must worry about two things: (1) the allocation of storage 
% for arrays, and (2) calling the appropriate TRAID routines to perform the 


desired operations. 
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2.0 TRAID PRIMER 

This section is designed to acquaint the reader with the basic con- 
cepts of TRAID. It shows simple solutions to simple problems; a more 
complete description of TRAID’s capabilities is presented in the next sec- 
tion. 


2.1 FUNDAMENTALS 

To illustrate the most fundamental aspects of the use of TRAID, let's 
choose a typical-but-easy problem and develop a complete solution for it. 
(Note that this problem is too easy to demonstrate the power of TRAID; it 


is intended only as a pedagogic device.) 


As a sample problem, say we are given the positions and velocities 
of an object and an interceptor, and we want a present-estimate of time- 


to-closest-approach. Expressed in vector equations, we have 


a 


re |(p3 + vor) - (Pi + Vit) |? = 0 


2(P6 - pt)-(vd - vi) + 2(V0 ~ vi)?t = 0 


pw _ 26 = PT) + (WO ~ Vi) 
(vd - vt)? 
where t is the time (from now) to closest approach and the (current) values 


of position and velocity for object and interceptor are Pd, PI, vo, and VT. 


Now, to translate from an algebraic expression to a FORTRAN program, 
we need to do three things: 
define position and velocity vectors 
calculate the vector differences 


calculate the dot products 


The first of these is ordinary FORTRAN lore. Let's declare 3- 
component vectors PO, PI, VO, and VI, and for the differences, DP and DV. 


The corresponding DIMENSION statement is shown below. 


CT Tr wy 


As to calculating the vector differences, TRAID offers a subroutine 
for just this purpose. Its name is SUBVEC, and it takes three arguments, 
each a 3-component vector -- SUBVEC subtracts the second vector from the 


first, and stores the difference in the third. 


The easiest way to calculate dot-products is to use the DOT function, 
which cakes two vectors as arguments and returns their dot-product as the 


function value. 


The resulting code is: 
DIMENSION PO(3), PI(3), VO(3), VI(3), DP(3), DV(3) 
CALL SUBVEC (PO, PI, DP) 
CALL SUBVEC (VO, VI, DV) 
TIMTOCA = - DOT(DP,DV)/DOT(DV,DV) 


So much for the calculations. We have yet to arrange for reading-in 
the position-velocity data and for printing TIMTOCA. TRAID provides sub- 
routines for reading and printing data, and it is recommended that you use 


them exclusively. 


The simplest input routine is IN1; it reads data from a standard- 
format card, one value per card. You CALL IN] (X,N), and IN] will read 
N cards, storing successive values in (the list) X. A more convenient 
routine for vector-data, however, is IN3, which reads three values from 
each data~card. These cards may have any descriptive text punched in 
columns 1~40, and should have three values punched in col. 41-50, 51-60, 
and 61-70 in F10.0* format. 


A simple printout routine is RITEF; its use is: 
CALL RITEF (10HNAME/DESCR, VALUES, N) 


When called in this fashion RITEF will print the 10-character 
NAME/DESCR, followed by N VALUES, in F-format. If the numbers in VALUES 
are integers, you should call RITEI, which prints I-format. 


If you punch the decimal, of course, you may place the value anywhere 
in the field. 


A a ET RI stg sent tr 


A more flexible printout routine is subroutine WRITIT. When using 
WRITIT, you compose any string of BCD characters (like an ordinary sentence) 
which contains short strings of zeros where you want successive VALUEs to 
be inserted. Then you CALL WRITIT(VALUEL, INDENT, nHBCDSTRING-), where 
INDENT controls the left-to-right position of the printed line. You con- 
tinue to CALL WRITIT (VALUE 2) and CALL WRITIT (VALUE 3), on until you have 
provided all the VALUEs needed to fill the zero-strings -~ at this time the 
line is printed. An example of WRITIT's use appears in Sec. 2.5. 


Two subroutines are provided for skipping lines and pages on the 
printout. They are CALL LSKIP(N) to skip N lines and CALL HEAD (N) to 


write a page~heading and skip N lines on a new page. 


There is one more requirement -- and this is important -- all programs 
which use TRAID must make a special call of HEAD at the beginning. This is 
of the form: CALL HEAD (nHPROGRMNAME PROGRAMDESCRIPTION-). When called in 
this way, HEAD reads a data-card and performs several "supervisory" tasks: 


o Establishes time and page limits (to be enforced on later 


calls to HEAD) 

Establishes physical constants and unit-scaling factors 
Establishes control parameters for the computation -- 
e.g., choice of integration methods 

Prints a title page, including the name and description 
of the program, and the number and description of the 
problem now being solved 

Assembles a page-heading, to be printed at the top of 


each page of printout. 


Shown on the next five pages are: 

A program to solve our closest-approach problem, complete 
with input and output and special initial call to HEAD, 
and arranged to allow multiple cases to be computed. (The 
call of PREDATA causes the images of the data-cards to 

be printed.) 

. A typical data-deck setup 

i The resulting printout. 
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000003 
600003 
000004 
000006 
000010 
noon? 
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000024 
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000042 
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22 STATE VECTORS 


Trajectory computations in TRAID make use of (augmented) state 


vectors -- 10-component vector defined as follows: 
STATE (1) = time in seconds 
STATE (2-4) = position (x-y-z) in neteee: 
; STATE (5-7) = velocity (x-y-z) in severe feacona 


oe: isa nbs * 
STATE (8-10) = acceleration (x-y-z) in meters feecona- 
STATE(1) is the time when the body has the given position, velocity, and 


acceleration. 


These are measurea in a coordinate system which is determined by the 
control parameter KOORD -~ see the figure on the next page. For anv 
TRAIN computation that depends on the altitude, KOCRD must be set to the 
appropriate value: 0, 1, or 2. Normally, KOORD is set to zero (by the 
special call of HEAD); to change it to some other value XN, you simply 


CALL SETKORD (N). 


There are TRAID routines designed especially for the inrut and out- 
put of state vectors. They allow a variety of (external) coordinate sys- 
tems; they allow the use of BCD names identifying the vectors; and they 


will handle all or specified parts of the state vectors. 


The state-vector input routine is STIN; it does the following: 
: Reads data-cards, extracting any or all of: 
BCD names of vectors 


time-component 


: position-components 
velocity-components 

Transforms coordinates if necessary 

v : Stores the vector in an input array 


Prints the contents of the data-cards 


rn nr 


* 
Instead of meters, these will be in feet if KUNITS has been set, on the 
HEAD data card. 


ORI ee pee 


ee rn ae eee 


0) EIT EOIN TTD EN 


,OkO = 2 
GRAVIT origin at earth's surface 
flat eartr 


LS RUGHE 


erigin et eartr s Surface 
soherica. eart” 


GRAVITY 


NORTS POLE 


os 


— EQUATOR oS, 
ITLuE 


GREENWICH 


ZERO LONGI TUE 


KOORD = 2 
origin at center of earth 
spherical earth 


LONG 


The state-vector 


mates if necessarv and 


printout routine is STOUT--it transforms coordi- 


prints BCD names and/or time and/or position 


and/or velocity and/or acceleration 


The arguments of sil and STOLT are the same: 


e lliLe is a BCD title to be printed as indentification 


° KFORM is a control code (of 6 digits or less) which 


specifies the external coordinate system. For example, 


KFORM = Llu veans position plus velocity in x-y-z coordinates 


and x FORM = 


+UJ in@ans position only in alt-long-lat coordi- 
nates. 


° NAMES is the location of an array of BC: vector-names. 


lf the thousands' digit of KFORM is zero, a dummy argument 


mav De used tor NAMES. 


e STATES is the location of the state-vector array 
e NST is 


Lf NST 


the number of state vectors to be input or output. 

is negative, title and column-headings are not 

printed. Ii NST is zero, the title and column-headings 
are printed but no vectors are processed. 

In order to illustrate the use of state vectors in TRAID, let's 


look at a noisv-radar-observation problem. Let's provide, as input 


data, a radar location, a target-object position and velocity, and a "state 


vector of radar errers. Then we'll print a set ot ten noisy position- 


‘ velocity observations. jfhese observations are produced bv subroutine 


RADAR, as shown in the program listing on the next page. 
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ck te ee 


cennge amen s 


PiiMFA2 © AN PB? CEMUNSTRATION FOR PRIMER 2 © 


a waraR i CON « 


POSITIGN COORDINATES 


x Y Z 
aM xe aM 
0.000 502000 1.000 


OA YECT TRAS © 


PCSITION COORDINATES VELOCITY COOROINATES 
Y Y 


a 4 x z 
KK nM KM M/SEC msSEC M/SEC 
0-000 100.000 20-000 02000 =70000000 30000000 
RADAR ERWORS = 
POSITION COORDINATES VELOCITY COORDINATES 
BANSE AZIMUTH ELEVATION MAGNITUDE AZyMUTA FLE VATION 
xe DEG DEG M/SEC DEG OFG 
16000 12000 12000 2-000 3.000 3-00f 
WOLSY CBSERVATIONS © 
POSITION COORDINATES VELOCITY COOROTNATES 
x Yy r4 x Y 2 
xe KM KM M/SEC MsSEC M/SEC 
2596 99.108 200306 2630255 «— $ E97S652R 1020 THE 
2628 100,566 19.568 46629) ©7017 .492 907.306 
0307 99.928 1907346 41.178 #7003.607 #972 269% 
~ 092) 100.936 200148 1272756 ©7015 2995 #902 049K 
225258 99,896 202398 1132983 $6990.36) -105a.557 
wlLet4o 99.405 200392 1724298 ©6904. 509 ©)0872290 
“0928 99.957 20-541 702946 ©$6993.A65 -) 058.494 
2773 100.376 19-985 $2999 ©7008.917 962+ 38A 
123946 101.221 192380 $123.496 ©7029e517 8390772 
-103264 99.308 182406 291-621 7008628 900.660 
4 
8 
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2.3 TRAJECTORY INTEGRATION 


TRAID integrates trajectories with the following features: 


A SE GE mg ore oer 


ae thrust, mass, lift, drag, maneuver limits, and delayed 


response, for up to 3 stages 


: 252 guidance commands, either dynamically computed or pre- 
selected 
3. termination-rule defined by the user 


For example, vou might (1) define a vehicle with lift and drag only, 
(2) input some preselected maneuver commands, and (3) ask for the vehicle's 
trajectory to be integrated until it reaches some specified altitude, or 


some specified speed, or some specified heading, or some other criterion. 


Or you might (1) describe a rocket with three stages of thrust, lift, 
and drag, (2) provide a pruportional-navigation routine which computes 
guidance commands toward a target object, and (3) request integration of 


the trajectory for a specified time interval, or other criterion. 


In both of the above examples, a single call of one subroutine will 
produce the requested trajectory, as a collection of successive state vec- 
tors (the first vector represents the start of integration, and the last 


represents the state where the termination~condition is attained). 


The rocket description is contained in an array usually named C and 
dimensioned 3 (stages) by 16 (descriptors). These descriptors, again, are 
parameters such as thrust, mass, lift, drag, response, etc. fcr each stage. 
Another array, usually named CP and singly-dimensioned, contains staging 


schedule and guidance information. 


‘ There is a TRAID subroutine (FLIN) designed to read information from 


data-cards, and store it in the C and CP arrays. FLIN takes 4 arguments: 


TITLE, a BCD title, ending on a dash 


MODE, a three-digit number 


a. first digit controls reading of thrust/mass data 


b. second digit controls reading of drag/response data 
c. third digit controls reading of preselected 
meneuver commands 
: Cc \ rocket~descriptor arrays 


° CP where data is to be stored 


For data of types (a) and (b) above, a more complete description ap- 
pears under Subroutine LIN in Sec. IV. In the case of type (c) data, for 
preselected maneuvers, vou punch on cards the capitalized entities in this 
command: “pull ANC g's (in a direction normal tu velocity and described 
by PSI and PSIDOT) until parameter NAME reaches a value of VAL." Typical 


examples are: 


ANC PSI PSIDOT NAME VAL Interpretation 
25 +90° D) VEL.EL. -45° "Dive at 25 g's until 


you're heading at 45°," 


10 0 60 deg/s TIME INT 6.0 "Turn a 10-g barrel roll 
in 6 seconds." 


Any number of the maneuver commands may be stacked ~- they will be performed 


sequentially. 


The TRAID routine which carries out the trajectory integration is sub- 
‘ routine FLIGHT. When you call FLIGHT, you must provide: 
| . C and CP arrays to describe the rocket and its guidance, 
| and the STATES array to be filled with subsequent state 
e vectors. 
e the names of three other subprograms: 
° One to calculate commanded acceleration -~ e.g. FLIER 
° One to calculate actual acceleration -- e.g. FLAC or RV 


* One to interpret the termination condition -- e.g. FOAL 


integration control parameters: a time step, a parameter 
name and value for termination, and NOUT (which causes 

a vector to be inserted into the STATES array every NOUT-th 
time step). 

. NST, which serves two purposes: it advises FLIGHT how many 
vectors you -xpect will be stored in STATES, and it tells 


your program how many vectors actually were stored there. 


To illustrate the above points, the program on the next page calls FLIN to 
input a rocket description and again for a set of maneuver commands -- 
then it calls FLIGHT to integrate a trajectory and STOUT to print it. 


Notice these items: 


af The dimension of CP must be 20+5*N to allow N maneuvers. 

* FLAC, FLIER, and FOAL must appear in an EXTERNAL statement. 

® An initial position and velocity must be set up in STATES. 

: Before calling FLIGHT, NST is set to the size of STATES. 

. After returning from FLIGHT, NST contains the number of 
“rectors. 

. FLIGHT is asked to integrate in 0.l-second steps until the 


altitude reaches 40 km, and to insert every 10th vector in 


the STATES array. 
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2.4 ORBITS 


TRAID offers an alternate representation of position-velocity -- ; 


) namely Keplerian erbits. It is true that most of the geometric calcula- 
tions and all ot the trajectory integration take place in the rectangular i 
E § 
« State-vector coordinites. bven so, there are some problens and parts of f 


problems which can be scived sore easily by using orbits; and in many 
problems the easiest solution involves both state vectors and orbits, and 


) frequent interchanges between the two. 


TRATD uses an orbital-element set to describe a Keplerian orbit. 
te Sv lonvention, this is a l0-component vector; it can describe either an 


ellipticai or hyperbolic orbit. 


The subroutine which reads and/or prints orbital-element vectors 
is subroutine ORIN, It takes 5 arguments: 


i . A BCD title 


A 3-digit control code. If the first digit = 0, ORIO prints 


but does not read data-cards. If it = 1, ORIO does both. 


MRE wea, 


If the second and third digits = 0, the 3rd argument may be 
a dummy. 

: A list of names (may be dummy) 

Location where orbital elements are steered 


Number ot orbital-element sets to be input/output 


The typical uses of orbits fall into two categories: transforming 
z from state vector(s) into orbit, and transforming back from an orbit to a 


Pa state-vector. 


For the first of these transformations, TRAID offers two options. 
Subroutine ORBL creates an orbital-element vector from the (1) time, 


(2) position, and (3) velocity in a state vector. 


The second option is ORB2, which passes an orbit through the positions 


| 
i of two state vectors (using the time at one of them) and satisfies one other 
if 


requirement: 


-- 


Flight-time is a specified value 


or Flight-time excess above minimum-energy is a specified value 

or Speed at one of the positions is a specified value 

or Velovity-elevation at cone of the positions is a specified ‘ 
valve 


To "reenter' from orbit into state-vector form, you call ORBP with one 


extra parameter. You mav choose between: i 


true anomaly (MODE = 0) 
radius (increasing) (MODE = 1) 
radius (decreasing) (MODE = 2) 
time (MODE = 3) 
altitude (increasing) (MODE = 4) 
altitude (decreasing) (MODE = 5) 


If you want only the time when one of the above parameters reaches a given 


value, you may use the function ORBTIME. 


Some illustrative uses of the orbit routines are shown in the example 


in the next section. 


2.5 EXAMPLE 

On the following pages is a sample program which exercises most of 
the features discussed throughout Secs. 2.1-2.4. This program finds the 
closest approach between an ICBM-type object and an interceptor which is 
launched when the object descends to some specified altitude, flies a pre- 
selected maneuver until it reaches 50 km altitude, and then goes "into 
orbit.” During the search for closest approach, both bodies are described 
by orbital elements -- these must be converted to state vectors for each 


trial solution. 


Notice the following items in the sample program: 
- There are many 10-vectors; ST is dimensioned 10 x 20, pro- 
viding room for 20 successive states; CP is dimensioned 45, 
allowing 5 maneuver phases. 


‘ FLAC, FLIER, FOAL must be declared in an EXTERNAL statement. 
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oe 


. The use of Keplerian orbits, in conjunction with atmospheric 
; flighe, dersands that KOORD be set to 2. 
The object's trajectory is to be determined by launch and 


impact-points, and a MODE-VALUE pair for ORB2. 


Wit ada Ng A 


The time for interceptor launch is determined by an input 


Bie treet oat coe cece * 


a) object altitude. The object's orbit time is reset so that 
this altitude is reached at time = zero. 

: e ; The interceptor's performance, its preselected maneuver, and 
: . its launch site are read in, and its clock ST (1,1) is set 

to the launch time = zero. 

‘ Before calling FLIGHT, MAX is set to 20, the size of the ST 

array; FLIGHT inserts every 50th vector into ST and stops at 

50 km altitude with MAX now containing the actual number of 

vectors in S1. The trajectory is printed and a set of orbital 

' elements is produced. 

In searching for closest approach, we start at the time the 


@ interceptor reached 50 km and we will search until we find 


Pee WE At 


positions within 1 second (TCLOSE) of closest approach. 

For each iteration, we obtain state vectors (55 and Sl) from 
each orbit, compute the differences DP and DV, and find the 
current TIMTOCA. If this is greater than TCLOSE, we reset 
TIMTRY and repeat. 


nee 


Soares 
e 


When we find a TIMTOCA which is less than TCLOSE, we advance 
the positions linearly for that time, find the position dif- 


“ ference, and print the time and miss~distance. 


There are three routines (which appear after statement 2C) 
which have not been described yet: 
®* CALL VECLIN (a,A,b,B,C) produces a 3-vector C which 
is a linear combination of the 3-vectors A and B -- 
l.e., C= aA + bB 
i * X = XMAG (A) returns the length of 3-vector A. 
¥ * CALL CEASE (12HBCDMESSAGE-) prints a BCD message, 
@ then prints a run-ending page, and stops the program. 
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f ® 3.0 TRAID GENERALITIES 


The remarks in this section are intended to 
Sec. 2.0Primer, which was deliberately simplified 
Sec. 4.0Subroutine Listings, which do not present 
To state it another way, if Sec. 4.0 describes the 


making up the TRAID family, this section provides 


fill the gap between the 
and incomplete, and the 
a coherent view of TRAID,. 
building block routines 


the blueprint for assembl- 


ala oatd 


ing them and the mortar which holds them together. 


. On the following pages, arranged by subject matter, are described 
. TRAID's operating principles, some comparative evaluations, and hints for 


use. 


oe ore Dan 


| 3.1 PROGRAM CONTROL 
t Operating principles are: 

| : Execution-time-limit and printout-page-limit may be input via 
data-card; exceeding either limit causes run termination and 
(normal) exit to the monitor. 

Endfile marks on the input cause run termination and (normal) 


exit to the monitor -- this feature may be overridden if you 


nat haat 


wish to continue computing after an EOF. (See function 


CHEKFIL.) d 


q ‘ Error conditions arising during TRAID calculations -- such as : 
| 7 2 . 2 4 4 
i” failure to converge -- cause run termination and (normal) exit 3 
} i 
i to the monitor; this feature may be overridden (if, for example, 3 


you wish to abandon the current case and proceed to read data ; 


} for the next case) -- see function MISTAKE. 
° 
COMPKG. This is a main program which reads FORTRAN cards from the " 
input file, inserts COMMON packages among them as needed, and writes the ; 
. FORTRAN-plus-COMMON on a file named MERGED. See Sec. 5.1 for typica! : 
’ deck setup using a COMPKG ‘control’ card. 


PREDATA/MORDATA. 1f you want to have all your data cards printed at 
the beginning of your run, call PREDATA. When using PREDATA, you must call 
it before any other TRAID input routine, and you must call it only once. 
t 4 Normally, your main program card will specify TAPF5=INPUT; if using : 


PREDATA, this should be simply TAPES. 


teeta es cool : ae a alaiininibiidalishits ites cnvicit saccansaas meinen 


lf you want to be able to edit your data-cards (in order to run suc- 
cesive cases), you may call M@RDATA. This routine writes an edited set of 
data-cards on TAPE5; editing may be accomplished by matching columns 1-20 
of an old data card, or by more explicit requests to insert or delete spec- 
5 ified numbers of cards. When using M@RDATA, your main program must provide 
a buffer for TAPE4. fie complete rules for this data~editing routine may 


be found under PREDATA/MB@RDATA in Sec. 4.0. 


HEAD. The first TRAID routine called -~ except for PREDATA -- 
should be CALL HEAD (nHPKOGRMNAME PROGRAMDESCRIPTION-). At this call, 


HEAD reads a data card containing: 


e A run number and run description, which together with PROGRMNAME, 
PROGRAMDESCRIPTION, the time, the date, and the page number, 
are written at the top of each page. 

° Parameters which control the calculations--i.e., KSTINT and 
KUNITS (see Sec. 5.6). 

> ° Optional execution-time-limit and printout-page-limit--these limits 

are enforced by HEAD when HEAD is called upon to skip a page 
on the output. If no value (zero or blank) is input for (either 
ot) these limits, then HEAD will not cut your job off; this does 


not affect the monitor system's cutoffs. 


It is on the occasion of this special call that HEAD sets up values for 
the physical constants used in the calculations, prints a title page, and 


sets KOORD to zero (see SETKORD, next). a | 


SETKORD. TRAID's interpretation of coordinate systems depends on 
KOORD (see the figure in Sec. 2.2). There are in TRAID three areas of 
sensitivity to KOORD: 

: Atmosphere models, embodied in DNSITY and GRAV 

. Calculation of altitude, as in ALTF 

7 Interpretation of radar measurements in subroutine RADAR 
Since the special call of HEAD causes KOORD to be set to zero, you have 


4 to CALL SETKORD (N) to give it a nonzero value. 


CLAN a Nhmamniintiyanbbicwnce + 


SETK@RD also offers entry points convenient for (re)setting KSTINT, 
KUNITS, AND WBODY. 


CHEKFIL. All TRALD input routines use function CHEKFIL to test 
@ for EOF on input. Normally, CHEKFIL will exit upon an EOF. To disable 
this automatic exit, CALL SETIFEF (1); to reenable, CALL SETIFEF (0). 


te 


While the automatic exit is disabled, you ought to check for the occur- 
rence of an EOF by looking at the number of endfiles, N = NUMFILS (DUMMY), 


after every call of an input routine. 


WHEN, This routine prints elapsed computer time. 


MISTAKE, All TRAID routines with error conditions call TRADERR, 
. which prints an appropriate message and exits to the monitor, To disable 
this automatic exit, CALL SETIFER(1); to reenable it, CALL SETIFER(0). 
While the automatic exit is disabled, you ought to check for errors by 
consulting the number of errors, N = MISTAKE (DUMMY). [You might use this 
feature if you wish to proceed to a subsequent case even if there is an 


error in a previous case. ] 


Q8ERROR. This routine is called by TRADERR to write the trace of 


calling programs; it prints the name of, and relative address in, each 


> subprogram in the chain back to the main program. 


CEASE, Subroutine CEASE provides a convenient way to print a mes- 


Sage,a run-termination page, and exit to the monitor -- all in one CALL 


PUT mre cS 


statement, Example: 
IF (NIRIES .GT. LIMIT) CALL CEASE (19HbITERATION FAILEDb-) 


3.2 INPUT 
Operating principles are: ad 
° * When data is read in, it is normally printed immediately in 
its original form. Thus, in normal operation, the printout 
shows exactly what data was used, 
a by If an endfile is found, normally an exit is made to the 
’ monitor -- this may be overridden, however. 
Data-cards are typically divided into 10-column fields, with 
BCD description on the left and numbers on the right, frequently 


beginning in column 31 or 41. 
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tq . All numbers are punched with decimal points: they will be 


converted to integers internally if appropriate. 
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If 'OLD DATA' appears in column 1-8 on a data-card, printing 

and further reading are suppressed and no data is transferred 
into core (until the next call of an input routine). 

If ‘END DATA' appears in column 1-8 on a data-card, this stops 
open-end reading (such as the reading of already-named variables 


in UGEILT and the reading of preselected maneuvers in FLIN). 


Some input routines accept a BCD title (which must end on a 
dash) to be printed above the table of printed data. 

Some input routines will read and return names from the left 
ends of data~cards -- typically, the NAMES are read from 
column ]l-30, and they are stored in an array NAMES (5,LINES). 
Some input routines allow their printing to be suppressed, 
Some input and output routines use a SETUP array, which contains 
an output format and column-headings for up to ten 12=column 
print-fields. This SETUP array is assembled by subroutine 
OUTSET, and used by INCOL and OUTCOL., If you call OUTSET, 
you must provide a SETUP which is dimensioned (82), 

All input is read from logical unit 5 -- except for PREDATA, 


which reads from INPUT and writes cn 5. 


Subroutine FLIN reads rocket-performance cata and preselected 
Its operation is controlled by the 3-digit MODE, as follows: 
The tirst digit triggers the reading of either 0, i, 5, or 7 
cards describiny thrust, mass, and burning schedule, 

The second digit, for values 0 through 4, triggers the reading 
of either 0, 2, 5, 9, or 11 data-cards describing axial and 
normal force coefficients, maneuver limits, and response. If 
this digit is a 5, FLIN reads 11 data-cards (as for 4), then 
calls IN] to read an integer N from another data-card, and 
proceeds to read N more data~cards in the standard FLIN format. 

If this option is used, the C array should be dimensioned (3, 
16+). 

The third dipit triggers the open-end reading of preselected 


Maneuvers. Reading continues until (1) a mark is found in 


; e column 71-72 of the last card of actual data, or (2) ‘END DATA' 
: is found in column 1-8 of the card after the last card of 
‘ actual data. To allow room for N phases of preselected man- 
: / euvér, CP must be dimensioned (20 + 5N): if there are no 
maneuvers, the dimension may be only (13). Note that it is 
a no longer necessary that the last maneuver have an unreach- 
able goal. 
a ° If you call FLIN for a rocket to be launched at a time other than zero, 
£ ° you should set CP(7) = desired Launch time, before calling FLIN. 
i 
j VYRLO. Subroutine URI reads orbital-element sets, which may be 
punched by hand (but rarely are), or which ORIO will punch for you if you 
call it with MODE = 2Xx, 
STIN., The TRAID routine which is designed to read state vectors is 
subroutine STIN, It will read time/position/velocity data in a variety 
Zz of (external) coordinates, such as: 
x- y- 2 
r-68-96 
: R-A-E 
alt - long - lat 
The choice of coordinates is implied by the 10's and 100's digits (KP and 
KV) in the second argument of the CALL statement -- but these may be over- 
ridden by KP and KV punched in columns 71-72 of the data-card. States read ” 
e by STIN are always transformed to standard state-vector form--note that the 


converse of this rule does not apply to STOUT. STIN stops reading when it 
finds an "END DATA' card, and resets LINES to the actual number read. 

a4 In], IN3. Subroutine IN] reads N values, one value per card, and 
scales it by either (1) multiplier/divisor punched in column 51-76, or 


(2) unit-name punched in column 51-53. If the second argument N is posi- } 


tive, data is left in floating-point form; if negative, it is converted 
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to integer form and an asterisk is printed on the printout; if N is 

zere, one card is read and the BCD characters in column 41-50 are returned. 

Subroutine IN3 reads N values, three values per card, with no scaling ( 
Capability. i: the resuire.d number of values is satisfied in the ‘middle’ 


Of a card, the values or the rese€ of that card are ignored. 


INMV. Subroutine (MMV reads from one data card a mode and a velue 
--e.p., for use by ORF2 or ORBP. The value may be input in such units 


as DEG, NMI, ete., if the appropriate name is punched on the data card. 


INDU, This is « routine which sets up and controls ‘looping’ vari- 


ables, aA call ot INDO causes it to read data-cards which define the loop- 


a ee 


: variables -- the cards say "let X vary from A to B in a manner M," The 
manner M may be either linear or geometric progression or random; and the 

z loops may be nested or concurrent, Later calls cause INDO to advance the 
variables to their next values. Note that (1) if all your loops are random, 
INDO cannot terminate and you must stop it, and (2) the loop variables are 


real (floating-point) variables and they must be stored contiguously in 


core. 


UNAMLT/UGETIT. TRAIL provides a routine for the purpose of 


changing the values of variables named on data-cards ~- this allows you to 


x delay the decisioa on variables-to-be-updated until data-punching time. 


5. 
{3 


Subroutine UNAMIT has the following features: 
™ The ‘changeable’ variables must appear individually in calls 

to UNAMIT -- on these calls, UNAMIT saves the address of the 

variable and the name which appears in column 1-]0 of the data- 

a card, and of course, stores the value in the appropriate loca- 

tion. 

The changeable variables need not be contiguous, and the names 

on the data-cards need not match names which are built into the 

program -- they only need to be self-consistent. 

Data-values may be scaled by multiplier/divisor or by unit- 

name (as in IN1). 


Data-values ure converted to integers if the name in column 


1-10 begins with a letter I through N. 
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After UNAMIT has built its table of changeable variables, you 
CALL UGETIT to cause data-cards to be read, their names matched, 
their values scaled and stored in the appropriate locations. 

zt : Either 'END DATA' or 'OLD DATA’ in column 1-8 on the card after 
the iast card of actual data, will cause UNAMIT/UGETIT to stop 
reading and return, 


Up to 30 changeable variables are allowed. 


INDEC/ININT/INALF/INNEW. TRAID has a format-free input routine which 
allows the mixing of descriptive text with numerical data, and also allows 


repetition factors. 


You call INDEC or ININT with an ARGument location and a quantity N; 
cards will be read and interpreted until N data-values have been accumu- 
lated and stored in ARG (they will have been converted to integers if you 
called ININT). Only legal number fields are considered; a number field 
will be ended on a blank or any illegal character, such as a letter # E, 
a second dot, etc. If a number field ends on a left parenthesis, that 
number is interpreted as a repetition factor, which applies to everything 
between that left parenthesis and the next right parenthesis that termi- 
nates a number field. On a subsequent call, INDEC/ININT resumes scanning 


: Pe the card-image where the scan was left off during the previous call, and 


reads new cards as needed to accumulate N values. 


v® 
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Entry INALF returns N words of BCD data from input cards. For 


each field, leading blanks are skipped, and the field ends on (1) the 
10th character, (2) a comma, (3) a blank, or (4) column 80 of the data- 


Card. { 


Entry INNEW sets a switch that prohibits the resumption of scanning 
of the current card -- so that on the next call INDEC, ININT, or INALF 


will read a new card instead of proceeding with the current one. 


TITLIN. This routine reads a title card in 8Al10 format; the title 
might then be used in calls to TITLER or SUBHEAD. 


INCOL., Subroutine INCOL reads data by a format which corresponds to 
the SETUP array assembled by subroutine OUTSET. Data may be converted by 
F, I, A, or O-conversion (see OUTSET). Note that INCOL assumes DIMENSION 
DATA (N, LINES) where N is implied by the contents of SETUP. If INCOL 
finds an END DATA card, it stops reading and resets LINES to the actual 


number of lines read. 


CHEKFIL. Function CHEKFIL checks for an endfile and (1) if there 
isn't any, returns, or (2) if there is one, then exits to the monitor or, 


if control flag IFEOF is set, returns. Typical usage is: 


READ (5,f) list ae 
IF (CHEKFIL (5) ) return control 
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OLDATA. This is the subroutine which TRAID input routines call, 


upon finding 'OLD DATA’ in column 1-8 on a data-card. 


3.3 OUTPUT 


of Principles of operation: 


& 


TRAID output routines count lines-printed and skip automatically 
to a new page when appropriate. The number of lines per page 
may be specified on the data-card read during the special call 
of HEAD (see Sec. 5.7); if not specified, HEAD sets it to 57 
lines per page. 

Some output routines accept a BCD title, to be printed, usually 
above a table of data. The first word of this title: 


If blank, causes no printing but a 2-line skip. 
. If 1H*, in some cases, prevents all the printing which 
would normally be done upon this CALL. 

If an ordinary string of characters, implies that this 

is a bona fide title which either is 1“ words long or 

is terminated by a dash (after the 10th character). 

Some output routines use the SETUP array which is assembled by 
subroutine OUTSET. The SETUP array contains an output-format 

and BCD column-headings. It is the mechanism by which TRAID 
produces its standardized tabular printout, which is diagrammed 
on the next page. By selective calls to OUTSET, any or all of 
the 10 data-fields may be used. 

Some routines allow suppression of all printing (see TITLE = 1H*, 
above) or suppression of parts of the printout, as follows: 

° TITLE = blank, described above. 

: LINES (the number of vectors to be printed), if negative, 
suppresses printing of TITLE and column-headings,. 

LINES = 0 suppresses printing of the data, but TITLE and 
column-headings will be printed. 

e All output is written on logical unit 6, except by SETPLOT 


which the calling program controls. - 


DIAGRAM OF TYPICAL TRAID 
ty TABULAR=PRINTOUT FORMAT 


120 print-columns 


e 
. BCD TITLE ENDING ON A DASH - 
cs HEADING HEADING HEADING 
ae FOR FOR FOR 
; FIELD 1 FIELD 2 FIELD 3 FIELD 10 
? 
xxx XeXX xx 
XXX Xe XX XX 
12 col> +12 col? +12 col+ 
NUMERIC CONVERSION IN EACH DATA-FIELD MAY BE: 
: E12.4 b+9,12344123 
rf? 
t E11.2,1X bb+9.12+123b 
«> F11.6,1X +999 .123456b 
<> 


F11.3,1X *#999999,123b 

F10.0, 2x +999999999bb - 
18, 4X +9999999bbbb 

I10, 2x +999999999bb 

012 777777777777 

2XA8, 2X bbABCDEFGHbb 

1XA8 ,A2,1X bABCDEFGHPQb 

(see subroutine OUTSET) 


} . Some routines allow vectors to be identified by 30-character 
7 names, from an array NAMES (5, LINES). 
6 - Routines are provided which you may call in order to skip 
lines or pages. 
: It is intended that in ordinary circumstances, your program 
need not contain its own WRITE statements and formats, but o 


should let TRAID do the printing. 


¢ 
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QRIO. Subroutine ORIO prints (or punches) orbital-element sets. 
N When printing, the vectors may be identified by (1) time in ORBEL(1), (2) 
» BCD name in ORBEL(1), or (3) the vector's sequence-number. When punching 
(KR = 2), ORIO produces cards which can be read-in when ORIO is called by 
o¢ a later program (with KR = 1). 


STOUT, STALE. Subroutine STOUT is designed to print a set of state 
vectors, either in standard form or a few variations thereof. The vectors 
. may be identified by 30-character names, and further identified by (1) time 

in STATE(1), (2) BCD name in STATE(1), or (3) the vector's sequence-number. 
Standard form state vectors may be printed in x-y-z coordinates (KP = 1), 


Printing of nonstandard forms is described in the following table: 


i Coords for Form of How STATE is transformed 
‘ Code Printout STATE from the standard form 
KP=2 r,0,¢ polar CALL STREP (STNDRD,STATE,+N,X) ; 
KP=3 R,A,E polar w iT " ! sot 
«> KP=4 alt,long,lat polar e si gs " vee 
tr KP,KV, (as for (special) User's program places a position- 
& KA>5 KX-5) deviation in STATE(5-7) or a 
velocity-deviation in STATE(8-10). 
KV=5 6,049 polar CALL STREP (STNDRD,STATE ,+N ,X) 


STOUT performs half the work of printing; it calls STALE to calcu- 


late the scaled-for-output vectors. 
LSKIP. This is the line-skipping routine. 


e HEAD, When you CALL HEAD(N), this routine (1) checks against the 
‘ time and page limits and exits if appropriate, (2) skips to a new page 
and writes a page-heading, and (3) skips N lines on the new page. If N is 


e zero, HEAD prints a run-end page, and returns. 


SUBHEAD. Subroutine SUBHEAD is called by HEAD just after printing 
the heading on a new page--and SUBHEAD prints any subtitles that have 
been established. The way to request a subtitle is to CALL SUBHEAD (n, 
> BCDSUBTITLE) where n is the subtitle number (< 3). Subtitles may be 
c changed or deleted by other calls of SUBHEAD. 
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CEASE. Called with a BCD message, CEASE will write the message and a 


’ terminate your program's execution. 


TITLER. This routine will print a message; you specify how many 
oe columns it is to be indented, and you may either allow (+) or suppress (-) 


centering in the remaining columns. 


COUNOUT. There are two special situations in which COUNOUT is 


? useful. 


Before your program prints N lines of output (on logical unit 6), 
calling COUNOUT(N) will cause the line-count to be updated, and the old 


page to be ejected and a new one titled if necessary, 


If your program is designed to input-compute-print for an unknown 
number of successive cases, and you want the next case's input to appear 
ae on a new page, then you should CALL COUNOUT(-60) before recycling to the 
am input routine, To illustrate: 
input common data 
10 input data for a case 
compute and print for this case 
CALL COUNOUT (-60) 


GO TO 10 


RITEF/RITEL/RITEA/RITEO. Designed for debugging-output purposes, 
this routine prints a 10-character description plus a specified number of 


| 
i 


data-values, converted by E-or-F, I, A, or 0 conversion respectively. 

This routine may be turned off or turned back on by calling RITEOFF or 
RITEON; and it counts and prints (1) the total number of times called 
and (2) the number of times called with the current 10-character name. 


OUT1, OUTN. 


or a dash) and either 1 or N data-values, you may call OUT1 or OUTN re- 


To write a short message (ending on the 40th character 


spectively, The output-conversion formats are chosen by the same code as 
used in subroutine OUTSET, In subroutine OUITN, sequential data-values are 
normally printed left-to-right, five values per line; if N is negative, 


however, they will be arranged top-to-bottom in five columns, % 
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RYDIT, WRITIT, These two routines can be used to insert and print 


a few data-values anywhere in a line of descriptive text. This text is a 

string of any BCD characters; the text-string should contain substrings of 

zeros where you want successive data-values to be printed. [These sub- 

strings may contain a decimal, but they must begin with a zero,] Both 

RYDIT and WRITIT proceed to encode a data-value into BCD characters, in- 

serting them into the text-string and taking proper account of the size 

of the data-value versus the size of the zero-string and the location of ‘ 
the decimal if any. One data-value is thus inserted into the text-string 

each time the routine is called; the text-plus-all-values is printed when 

it is completed. Differences in usage are: 

7 RYDIT requires that the mode (real or integer) of the data- 
value be specified; WRITIT deduces the mode from the apparent 
magnitude, 

WRITIT's interpretation of INDENT is the same as in TITLER; 
RYDIT does not center, but does simple indenting only. 

The text-string for RYDIT may be up to 120 characters long 

and it must end on an entire blank word; WRITIT allows 136 
characters in length and termination on either a blank word 

or a dash. 

RYDIT must be instructed "do not print, but wait until the 
text-string is filled" (by setting MODE negative); WRITIT 
searches the text for unfilled zero-strings and waits-or- 
prints automatically. 

RYDIT inserts the encoded data-values directly into the text- 
string, which is an array stored in the calling program -- this 
means that (1) RYDIT may be used to assemble more than one 
text-string concurrently, but (2) after a given text has been 
filled with data and printed, its zero-strings must be re- 
stored for subsequent use. WRITIT, on the other hand, saves 
the text-string internally; this means that it must not be 
used "in parallel," but must be allowed to completely fill and 


print one text before starting on the next one (or restarting 


on the same one). 
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WHEN. This routine will help if you are investigating the computer- ae 
time usage of a program. For both central processor and peripheral 
processors, it prints both the cumulative time (since beginning of job) 
and the elapsed time (since the last call). It also prints the sequen- 
tial number of this CALL, the name of, and the relative address in the 


calling routine. 


OUTCOL. Subroutine OUTCOL is the routine which prints arrays by the 
format, and with headings, in a SETUP array. OUTCOL assumes that the data 
is in an array dimensioned (N,LINES) where N is implied by the format in 
SETUP. 


OUTSET. This is the subroutine which assembles the 82-word SETUP 
array, containing an output format followed by printfield-headings for up 


to ten 12-column fields. 


3.4 VECTORS AND COORDINATES 
The TRAID routines which manipulate vectors are here divided into 


(a) those which work on 3-component vectors only, and (b) those which work 


on standard and nonstandard-form state vectors. 


1. 3-Vectors Only 
DOT. Dot-product. 
cae 


CROSS. Cross-product. 
CROS1, Cross-product, made unit-length,. 
UNITV. To make unit-length. 


XMAG. Length of vector. 


AXVEC, Sets up either Kth coordinate vector, or (if K'= 0) null 


vector. 
PROJ. Projection of a vector onto a given direction, 


SEPA, Angle between vectors. 
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SUBVEC, Difference between vectors. 
VECLIN. Linear combination of two vectors, stored in a third. 
VECSUM. Linear combination of two vectors, added to a third. 


AZF. Azimuth angle between vector and the x-axis. 


ELF, Elevation angle between vector and the xy-plane. 


LOCLAX, Sets up a triplet of unit vectors, aligned with given 


vectors -- e.g., local coordinates at radar on spherical earth. 


2. 10-Vectors 
ALTF. Altitude of a state-vector-position. 


GRAV. Returns gravitational acceleration (3-vector) for a given 


state-position. 


(RADAR. Produces a set of noisy radar observations, using as inputs: 
(1) a set of object positions, (2) a radar location, and (3) a radar~error 
vector. The error-vector is interpreted in radar polar coordinates. The 
radar observations are returned in either (1) absolute x-y-z2 coordinates, 
or radar polar coordinates, based on the radar, wich (2) azimuth measured 
from the x-axis and elevation measured from the xy-plane, or (3) azimuth 


measured from east toward north and elevation measured from the horizon. 


STREP, Transforms a set of vectors between standard state~vector 


form and polar-coordinate form (t, r, 6, ¢, r, 8, $, ¥, 8, $). 
SITEP,. Performs earth-rotation on a state vector, 


EULANG. Produces direction-cosines (10 by 3) from Euler angles, 


and vice versa. 


TRNSFM. Transforms a state vector, using (10 by 3) direction-cosine 


array. 
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3.5 TRAJECTORY INTEGRATION me 
| Principles of Operation: | 
. TRAID trajectory integration takes an initial state-vector 
and produces subsequent state vectors which are stored adja- ( 
s° cent to the first. 
° The user specifies the guidance commands -- these may be pre- 


selected or computed dynamically. 


: , The user specifies the vehicle's characteristics, such as 
thrust, mass, lift, drag, response to commands, etc, 
bs The user specifies a terminating condition for the flight -- 
e.g., "stop when ALTITUDE is 40 km," or "stop when TIME 
INTerval is 3.0 seconds." 
S The flight may be computed by rectangular (KSTINT = Q) or 
Runge-Kutta (KSTINT = 1) integration -- KSTINT may be input 
on HEAD's special data-card, or set by calling SETINTG. 
- = Flight may go forward or backward (DTS negative in call of 
Re FLIGHT) in time. (Each state vector has the appropriate time 


stored in its first component, of course.) 


FLIGHT, This is the executive routine for the trajectory integration 
task. FLIGHT must be called with three subprogram names as arguments: 
, GUIDE, a subroutine which accepts as inputs CP and the current 
: state-vector, and returns a guidance command in CP(1-3), a 
time-step, and an iteration counter. For this purpose you 
may use FLIER (described below) or you may write your own 
routine. 

% AXEL, a subroutine which accepts a guidance command in CP(1-3), 
and stores an actual acceleration in components 8-10 of the 
current state vector. Typically, AXEL accounts for thrust, 
lift, drag, and response, For this purpose you may use either 
FLAC or RV (described below) or you may write your own routine. 

i GOALFN, a function which returns the numerical value, from the 

< current state vector, of the parameter defined by NGOAL, For 
c this purpose you may use FOAL (see below), or you may write 


your own function. = 


When passing subprogram names as arguments (FLIER, FOAL, etc.) you must 


not forget to declare them in an EXTERNAL statement, 


In the call of FLIGHT, your program supplies a STATES array, (1) of 
which the first state vector is the beginning state of the trajectory 
integration, and (2) into which FLIGHT inserts (some of) the successive 
state vectors as the integration proceeds, The number of vectors inserted 
into STATES is controlled by DTS, NOUT, and of course the nature of the 
vehicle, the trajectory, and the terminating goal. 


Since FLIGHT integrates until your goal is reached, there exists the 
danger of over-filling the STATES array. To help diagnose this condition, 
it is recommended that before calling FLIGHT, you set the last argument 
NSTATE = the number of vectors that will fit in the STATES array. Then if 
the number of vectors inserted into STATES exceeds that initial value of 
NSTATE, FLIGHT prints a message -- notice that FLIGHT keeps on integrating 
until the trajectory~goal is reached (or until something gets clobbered by 
the overflowing STATES array). 


While integrating, there will be discontinuities in the trajectory -- 
one upon reaching your terminating goal, and perhaps others due to sudden 
changes in the guidance command. On these occasions, FLIGHT integrates 
past the discontinuity and then tries to fly back to it. It is necessary 
therefore that the GOAL-criteria have a continuous first derivative in the 


region of the discontinuity. To illustrate, a GOAL of velocity elevation 


= 90° will never work (because no heading has a velocity elevation of 99 + €), 


and a GOAL of velocity elevation = 89° will work reliably only if the pro- 
duct of turning rate (deg/s) and time-step (s) is less than 1°, 


If NOUT is positive, the state vector at the above discontinuities, 
plus the states at stage separation and ignition times, will be stored 
in the STATES array (as well as the periodic states every NOUT time 
steps). For NOUT negative, all of these extraordinary states are omitted: 
and if NOUT is zero, only two vectors are returned: the initial and final 


states. 


A simplified flow-chart of FLIGHT's operation is on the next page. 
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SIMPLIFIED FLOW CHART FOR FLIGHT 


s* 
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if CP(13)#40, 
call GUIDE to 

compute guidance 
command 


call AXEL to 
compute actual 
acceleration 


call DISCON: 
have we reached 
FLIGHT's GOAL? 

(GOALFN ) 


store vector 
in STATES 


store vector in 
STATES if 
appropriate 


call STINT 
to integrate 
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FLIER. This routine is designed to perform the GUIDE function for 
FLIGHT. It interprets the preselected maneuver commands (which have been 
read into the CP array by FLIN, called with third-digit MP = 1) in traject- 
ory coordinates, and sets up in CP(1-3) the commanded acceleration in x~y-z 
coordinates. The CP array must be at least 13 words long, even if no ma- 
neuvers are specified -- if there are N maneuvers, the length of CP must 
be 20 + 5°N. FLIER allows sudden changes in the guidance command; there- 
fore it uses an iteration-control flag to communicate with FLIGHT, and it 
returns a time~step DIE designed to reach the guidance-discontinuity. Note 
that FLIER uses FOAL. If FLIER finishes a maneuver-phase and finds that 
the next phase is undefined, it sets CP(13) to zero so that FLIGHT will 
not call FLIER anymore. 


FLAC, RV. These routines are both intended as AXEL-routines in 
FLIGHT -- both set up an actual-acceleration which includes commanded ac- 
celeration and vehicle performance characteristics, Subroutine RV computes 
this actual acceleration for bodies with no thrust and no response, but 
only drag; when using RV, the C-variable does not need to be an array but 
may be a single number--namely the ballistic parameter in kg/m? if 
KUNITS=0, or in ibe/£t? if KUNITS=1, (Hint--you can input this (in 


1b£/£t7) via IN] if you punch 'PSF' in column 51-53 on the data~card.) 


Subroutine FLAC provides a much broader vehicle-model than RV does. 
It uses the C and CP arrays to determine thrust, mass, maneuver limits, 
and response, and calls DNSITY, GRAV and SONIC for models of atmospheric 
density, gravity vector, and sound speed, For aerodynamic coefficients, 
FLAC uses functions CAXIAL and CNORML (which in turn refer to the C-array 
to compute the axial and normal-force coefficients). CAXIAL and CNORML 
are provided in the TRAID library, but of course you may write your own, 
Note that if you do write your own, and if you want to use a larger C-array, 
you should: 


. Remember to dimensiou C appropriately in your main program 


Not tamper with any of the existing C-parameters 1-6, 11-12, 
15-16 ° 


77 NO 


Input any extra parameters by calling FLIN with second-digit 
LD = 5 (don't forget the INl~type card which specifies how 


many extra parameters follow). 


FOAL, This function is designed to perform the GOALFN task in FLIGHT. 
It interprets the termination criterion NGOAL, then computes the current 
numerical value of the corresponding parameter, from the state vector. 
Note that FOAL is used explicitly by FLIER, while FLIGHT uses whatever 
function is specified in the CALL statement -- which usually is FOAL again, 
but may be another function if you want to write one. Most of the values 


FOAL returns are continuous functions; two kinds are not continuous: 


Value is aaa same | NGOAL criterion 


PHI 

ELEVATIO 
LATITUDE 
VEL. EL. 
ACC. EL, 


CAXIAL, CNORML, These functions produce axial and normal-force coef- 
ficients from mach-number, attack-angle, the C-array, and stage number. 


[See note under FLAC (above) about writing your own versions of these. ] 


DISCON, This routine chooses a time-step to try to reach a trajectory 


discontinuity. 


STINT. This is the routine that performs the integration by either 
rectangular or Runge-Kutta methods (controlled by KSTINT which is in common- 
block BASCON and is read from HEAD's special data-card). STINT does all its 
communication with FLIGHT through COMMON variables, not through arguments. 


DNSITY, SONIC, GRAV. These routines provide FLIGHT with information 


about the real world. Note that DNSITY and GRAV are sensitive to KOORD, 


ria 
as 


oe eee ee 


UPSTATE. This routine advances a state vector to a given time, 


using constant acceleration. 


TRPSTA. This subroutine interpolates between a pair of state 


vectors, to establish another state vector at a given time. 


3.6 ORBITS 

TRAID computes Keplerian orbits using an orbital-element 10-vector; 
the first of these is not used by TRAID except it may be input and output, 
the tenth is not used either, and the second through ninth form a redund- 
ant set describing the orbital ellipse or hyperbola. Routines which estab- 
lish orbital-element vectors (ORB1, ORB2, RORB2, and ORIO) always set up 


the whole redundant set, so that you may use any convenient combinaticn of 


them. 


ORB1, ORB2, RORB2. These routines calculate orbital elements from 
state vectors: ORBl uses STATE (1-7); ORB2 uses STATE1] (1-4), STATE2 (2-4), 
and one other parameter (for which you have four options); RORB2 uses the 
same inputs as ORB2, but you have six options. Both ORB2 and RORB2 store 
the appropriate time in STATE2 (1). 


ORBP, ORBTIME, ORBTIM, ORBPD. These routines use orbital elements 
as inputs and produce various answers in state-vector space. ORBP, with 
one other parameter (your choice of six), produces a STATE(1-7); ORBTIME, 
with another parameter (your choice of six) produces clock time; ORBTIM, 
with another parameter (choice of four), produces time-since-perigee; 
ORBPD, with time, produces partial derivatives (dimensioned 66) of STATE 


(2-7) with respect to ORBEL (2-7). 


3.7 MISCELLANEOUS 
DATEF, Returns BCD date in the form 'MM/DD/YY'. 


DTIMEF, Returns BCD time-of-day in the form 'HH-MM°SS', 


KALLER. Returns the name of, and relative address in, the N«th 


level calling routine. 


LOC2. Returns (address of 2nd arg) - (address of lst arg) +1. 
Typical usage: ‘ 


N = LOC2 (START, END) ‘ 
” CALL XMIT (-N, 0, START) 


LOCFL. Returns (field length) - (address of arg) +1. 


MLXER, XMIXER. Provide a convenient way to store integers in 


floating-point arrays, and vice versa. 


FDIV. Replaces an illegal operand (infinite or indefinite) with 


PACBIT. Packs and unpacks data. 


KSHFT, SHFT. Shifts a word to right or left, filling with zeros. 


INDMAX, INDMIN, INDEQU. May be used to scan a list, returning 
the index of: 


© INDMAX--the numerically largest element (the first one 
if there are more than one of the largest value); 
r INDMIN--the numerically smallest element (the first one 


if there are more than one of the smallest value); 


e e INDEQU--the first element that matches a given match word. 


RNV. Returns a quasi-normally distributed variable. 


KONVERG. This routine does the tedious part of an iteration, 
choosing new values for the independent variable, and testing for 
divergence and excessive tries. KONVERG is initialized by a call of 
" KONVSET; KONVERG stores its information in the calling program and 
C therefore you can stack any number of iteration loops. | 
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MINMIZE. This function works just like KONVERG except that it 
searches for a minimum value instead of a zero value, and that it 
assumes the function is parabolic (near the minimum) instead of linear 


(near zero). 


TFLITE, TFLYTE. Returns a flight-time between two points: TFLITE 
assumes that (1) thrust and drag may be adequately represented by an 
initial speed and time-to-attain-that-speed, and (2) gravity is constant; 
TLFYTE dispenses with the second assumption, but is only one-quarter as 


fast in execution. 


SONIC, DNSITY. Provide models for sound-speed and atmospheric 


density. 


TRPLATE. Performs N-point interpolation in a monotonic list. 


SETSCAL. Computes appropriate endpoints for the axes of a plot. 


SETPLOT, PLOTT. Produce plots on the printer. The primary differ- 
ences are: 
SETPLOT must be called with each x-y pair (data need not be 
stored in a table) and a character to be plotted (63 choices 
including blank). Characters may be added below or to the 
right of the previous character, which allows building up 
labels, Any number of traces may be plotted. SETPLOT uses 
a COMMON block of 663 cells. 
PLOTT works from a table of data, and assigns to each trace 
(up to 6 traces) a preset character. Once your data is 
arranged into tabular form, you can easily plot various 
columns versus each other, and you can easily plot a histogram 


instead of a single trace. 
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MATDIAG. Sets up a square matrix with a given vector on the diagonal. 
MATSCAL. Multiplies a matrix by a scalar. 

MATFLIP. Transposes a square matrix (into the same location). 
MATRANS. Transposes a matrix, storing in another matrix. 

MATINVRT. Inverts a matrix 

MATADD, MATSUB, Add and subtract matrices. 

MATMULT. Multiplies two matrices. 

MATIGEN. Calculates eigenvalues and eigenvectors. 


MATPANP, Produces a submatrix by selecting and reordering specified 


rows and columns from the original matrix. 


JACOBI. Is a more general form of MATIGEN. 


eee 


oo | : 


rental roe agg 8 Rin reo en marae eas. 
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4.0 TRAID SPECIFICS 
In this section are the complete FORTRAN listings of all TRAID rou- 


tines. Information which was too subtle for the Primer (Sec. 2.0) or was 
too detailed for the Generalities (Sec. 3.0), will be found in this section. 
? With luck, specific questions may be answered by the explanatory comments 
t 


in each routine; as a last resort the answers can be found in the source 


code itself. 


a> 


.? 
* 

*. 000003 

eecees 

000003 

000003 


FUNCTION OLTF (STATE) 


c SOURCE DATE 0720726 BRAND NEW CODE 
c 
c RETURNS ALTITUDE OF POSITION IN STATE (206) 
c 
CCOMPKGsBASCON 
COMMON /8aSCON/ 
I PROGRM, KPAGE? LINEs TOF DAY» RUN» 
2 MSGe FLAGe OaTEee MAXTMs MAXPG, 
3 KLASS» REGROUP. KUNITS, KSTINTs KOOROs 
@ eNOFILE 
c 
CCOMPKG se CONCON 
COMMON /CONCON/ 
1 Ple SRO» SLVe SMF, SKPe 
2 GacCe GCONe wpoor. AHOZROs TaOPT 


OLMENSION STATE (10) oP (3) 


CALL AMIT(3oSTATE (2) oP) 
IF (ROORV?]) Le2e3d 
1 a ® Pt3) 
60 TO 4 
2 P(3) @ P13) e¢RBOOY 
3 A & MAGI) -RBODY 
oarly oa 
RETURN 
ENO ALTF 


ALTF 


RUNTN(6)¢ 
MARLNe 
FEOF 


RAODYs 
marry 


cacsie7ibubviaasah groveranshsndpensrcmaneateaginedaingend 


| 


e00004 
000004 
6 


SUBROUTINE AAVEC (Kev) 
SOURCE DATE 6700726 @RAND NE@ CODE 


RETURNS VECTOR Vv @ UNIT VECTOR IN KReDIRECTION 
OR VS NULL@VECTOR IF K © 6 


DIMENSION (3) 


Vt1) av (2) 0403) 800 
TECH eGT 00) VIRDSLe 
RETURN 

END = AxAVEC 


one 


ee 


een yew 


é CAXIAL 


{ « FUNCTION CAXIAL (VMACWeALFeCok eco? 68,0311 
c SOURCE DATE 68.1205 USE FOIV WHEN NENONINATAR MAY BE ZERO 
c @OURCE DATE 68,0311 ADD SeTH ARGUMENT 
¢ @OURCE DATE 68.0310 CHANGE ORAG PROFILE TO PORM CreZon/s(¥MaCheoy) 
¢ SOURCE OATE 67,0707 ADAPTED FROM FLAC OF 66.0001 
4 c 
4 € RETURNS AXTAL oPORCE COEFFICIENT, AT MACH NUMAER VMACH AND 
© ¢c aTTAck ANALE af, POR & VEHICLE PESCRIBED BY THE KotH STAGE 
s c ar THE C aRnRay 
¢ 
oeeei0 AIMENSION C(3,16) 
eoeo ie DATA (VL0065) 0 (VZ202,0)0 (VaeIe3e (VealOs? 
¢ 
0000)0 1H (¥MaCH,ATeV1) 60 To 2 
o0ee13 CARTAL eC (me?) $ RETURN 
000048 2 TF (VMACH.OTeVE) 60 TO 4 Z . 
oecesi CART ALEC (Ke 7) o (VMACHoVII/IV2VIIO(CINeBOCIKe7)} S RETURN 
© PalVOo¥3) (VIZ) OPO TY (CK ODOC (Ke FIFI C (Ket oC (KelO)?? 08,1205 
vaPOIv( (20¥4) 7 (Z0V!2)) 66,1205 
WO (C (9B) OC IK 9G) OLVGCYIO(VZ0V} /(V3eV2) 
vol (xB) oPDIV IRs (V2eY)) 60,1305 
CART aLaZoPOly (Ks (vMACHaY)) 68,1205 
eetuan 
ww CART AL 


CEASE 


a? 
OUT TNE cease imesSage) 
Cc squact pate 6720726 grano NEF coo€ 
¢ 
C provioes & wessace tas 16 T1TLenm ANO an Ext? 
©. ¢ ait You nave tO po 18 CALL 1 
¢c 
Sad 99eo0s caLl TITLER pomESSAGED 
eqeceé call wean ie? 
eeueré catt exit 
eeeets end CEASE 


aru 


CHEKFIL 


FUNCTION CHEKFIL (LUN) 60.0209 
SOURCE DATE 68.0209 CONVERT TO COC weo0 
SOURCE OATE 66.1215 MOFILES1s AND IFEOFaO IN DaTA STMT 


CALLS EAST OR SETS EOF FLAG IF ENDFILE ON LOGICAL UNIT LUN 


eee NOTE = IF AN ENOFILE HAS BEEN FOUND, THEN CHEKFIC WILL 
EAtT UMLESS *1FEOF™ HAS BEEN SET NON@ZERO 


TO LEARN NUMBER OF ENDFILES WHICH WAVE BEEN REAOe 
WenUME TLS (putty) 

TO RESET TO | THE NUMBER OF ENOFILES, 
CALL SETNFIL GL) 

70 RESET TO & THE ENOFILE@CONTROL SWITCH IFEOF. 
Cale SETIFER IL) 


aannaananaanane 


8 
i 
3 
4 
$ 


6eeees COMMON /BASCONS 
1 PROGAM, = KPages LINEs TOF DAYs RUN, AUNID (6) 9 
2 MSOs FLAOs DATEs MARTMs HAXOG, MAXLNy 
3 KLASSe KOROUP, KUNITS.  KSTINTs  KOORD, wrEOF 
4 


‘ 600003 OATA (MOF ILE@O) » (IFEOR Gg? 6601215 
a 00000) EQUIVALENCE (AQ 3€0) 


c 
r CHERFIL O00 66.0209 
IFCEOFeLUN)D 30000 


10 IF(SFEOF) 30020 
20 CALL HEA0CO) 
Cale EATT 
30 NOF ILEONOF ILES) 6691215 
Getta 100 6920209 
40 RETURN 


aepenertce s o8 e200 

Ie € 

Stun Iu *s€0 68.0209 
Lal 


ENTRY SETIFE F 6A .0209 
FEOF CLUN 
RETURN 


ENTRY SETAFIL On20209 
NOF ILE ®LUN 
RETUAN 


‘ 900087 ENO CHEKFIL 


Meee 
a 
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FUNCTION CNORML (VMACHSALF eC oKeCP) 
SOURCE DATE 68-0311 A400 SeTH ARGUMENT 
SUURCE DATE 67.0707 ADAPTED FROM FLAC OF 66.060) 


RETURNS NORMAL ©FORCE COEFFICIENT, AT MACH NUMBER VNACH AND 
ATTACK ANGLE ALFo FOR A VEHICLE DESCRIBED BY THE KeTH STaGE 
OF THE C ARRAY 


OLMENSION C( 3016) 


CNORML@ALFO(C (Ko) 3) oALFOC {Ky 36)) 
RETURN 
ENO CNORML 


CNORML 


68.0311 
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COMPKG 
® 


Pe RY EERE AMI oa os 


i 


PHUOGL AM CUKEKG 


¢ SOURCE UATE 8906709 wKTTE ENROW MESSAGE AMIL MENGE CARL 1MAGES 
a? c SOukCE UATE @HelZ0S PRINT CCOMPKG CARD bUTh HEFORE AND AF IER : 
c Suuket Lott beeujul demas Sew COLL 
x4 
c 2SSIGN BUFFERS 
x TAPUTZaUu0Ue 
K NUTPLTB20226 
ry > ERGEIZ4000s 
> x TAPER) SOUIPUTs 
x TAPE] BINPUTs 
~ x TAPE 2BINPUT? 
:s x TAPE LUZMERGEC? 
g Cc 
+o e LOADEN VIA CONTROL Cane THIS PRUGRAM KEALS FORTHAN MACKUS 
aa e (6G = FACKAGES UF CUMMON STATEMENTS) FrUe SINPUTZ, TREN REALS 
¥ e FURTRAN SUURCE COLE FMOm 2INFUTZ»_ AND PRKTIES Teal SOUKCE COLE 
74 e CINCLLLING TRE aPPROPRIATE MACKOS) ON FILE =MERGEDSe 
es ° EaCM MACRO IS a SET OF CARDS PRECEDED RY A NAME“CARD ®ITH 
¢ e sMACKOMAME® IN COL lel0e AND FULLORED oY & CAHC wITH ZENDCOMS 
e IN CCL l=6e THE LaST SUCH ENDCUM CakD IS FOLLOWED BY A Caro 
e wITh etAGaAlLhs IN COL bebe AnD NEAT IS thE FORTRAN CODE IN wHtCr 
e TmE waChKCS ake TO BE IWSERTED. IN PROCESSING THIS CODE? 
e Tr— PRUGRAM LUOKS FOR CARDS WITH SCCOMPRGe® IN COL l=8e TAKES 
e TRE SMACRONAME® FwOm COL 9-180 ANDO WRITES THE CAWD IMAGES OF ThE 
° CORNE SPCALING MACHUs EMBEDDED AMUNG The FURTHAN CAKDS»s ON SMENGEDS, 
© UPON READING AN ENDOFHECUNDe Tht FRUGHAMN KEWINOS SMERGED® AND EXITS: 
° TYPICAL USE OF TRIS PROGRAM IS IMeLIED dy THE TYPICAL DECK SETUP ooo 
° CONTROL CARES JOD coe 
S e (REC 0) COMPKG, 
e HUNGS oo MERGES 
° LGO. 
° 1/8/9 
e PACHOS PLUS CUOMNAME 1 
IL . FORT MAA COMMON JAS AeYod 
E . (WEC 1) ENOCOM 
: e COMNaAME DO 
. COMMON /)/ relebek 
e INTEGER &K & EWUITVALENCE (Ooh) 
e *NDCOM 
° ENDALL 
4 ° SUBRUYTINE BU Eed) 
ry CCOMPRGsCOMNaME? 
e CCOMPRG se COMNaME I 
° OO 29 Lele10 
° ETCeee 
ry 7/8/9 
° Nava (EC 3) DASA CakDS 
* 
4 . 000003 : DIMELSTUN SRRAY( 4000) sLOM(E2) cL (O02) eNAME (H2) eRARU( 9) 
000603 CUMFON ARRAY 
‘ 000003 INTEGER Amhay 
000043 DATA (MAAARAY &® 4c Ou) 
; 0u0nI3 LIME STOR mESG (O07) 
bs 060093 NOTAIPESE s 69.9709 
# 1 Jib FanITyeCheCR OW CEB tkrUT ’ 64.0704 
4 ss 2 31h COMMON STATEMENTS EACEEL AkNAYs 69.0709 
X 4 3 JIM BAWITYeCHECK UN SOUNCE INFUT 4 69.0709 
¥ Y 31M AtVEH MEARD OF (Name) ’ 69.0709 
i 5 31H CUMPRG = (NAME) INSERTED oe 69.0709 
% 6 31m PARITY*CHECK UN COMMUN INPUT 6 6900709 
| 7? 3ln TOO MANY COMMON PACKAGES N 69.0709 
} c 
‘ 600003 Luve}o § AUNS®] § BUNS? 
eo000ue Kae] 
000007 \ooM s 0 
8 000019 ARCLEBH) CUFPRG 3S BRITE (6le)) KECU 
000017 CALL APITCLOOATEF (KDATE) oRUATED 69.0709 
c 
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000024 an NCOMaACUMS) 
000026 #608) ‘ 
000927 TE CNCUMeGT«O2) GC TO ble 
000032 LUM (ACCOM) BRACB 
000036 Ananae? f 
000035 PC2 REALINLNG 1) CARMATIAD (RB ROK) 
000050 1 FOkMaT(B410) : 
000050 TF (EMF eNLN) 8000106 
00053 166 GL TO 1106201004%e1 064) oK6Y 
000062 1062 IF (ARRAY(AA) EGeBRENPALL =} GU TU 110 ! 
o0006s 10oe KReKae? i 
000067 If (KR eGTemAABMAY? GO TC B02 
000072 TF (CARMAY (RA) cEGeGMENCCUM 7 GU TO lle 
000074 1])0 Kasnne) 
000076 ROsKu eg 
900077 R003? 
000100 GO To 1Uu2 
000101 Ll] LRE(NCOM) anben 
000106 LELOw (NCC) 98 
Q001u5 RAME (h.COm) BARMAY (LD 
090107 IFtNCOmet ted) GU Tu 116 
000111 DU lie sCs2encor 
000113 1F (NAME (GCob) eGo NAME (NCOMIS NAME (yCel BH) 
000117 Lhe CONTINUE 
000122 TLS Ik tnGGeEGe3!) 11e6e10 
c 
000127 110. NCOMeACOMN=] 
c 
Cc 
000131 30 MEAD (MUN, 2) KARL 
000137 2 FURMAT (ARS TALC #A2) 
000137 If (EOF oR) 60038 
aca1e2 Je wwJTE(LUNe2) KARO 
000150 IF (RaMCeNE CEMCCOMPKUS) GU TO 30 
000152 IF (NCOMeL Tes? GU TO 3h 
000155 C0 34 NE} eNCCH 
000156 TF (NAME CR) cEwemahu(2)) GO TU 00 
000140 3o CUNT TRLE 
Ou0142 3m MESGi30G) & RARL(2) 
O001t4 ehL[ Tr (hled) (mESGCINe@) ehalee) 
000176 wrth (LURe]) (MESO( Ree) Kel 0g) 69.0709 
0002u7 bu Tt 30 
C 
000210 oy Ka = LOwth) 
000212 we MESG(205) = NAMEN) 
o00216 wnITTE(6lel) (MESG(K95S) Red se) 
900226 68 Kb & Khe? 
000230 WRITE (LUN]}) (ARRAY (KR) oKanAent) 
0002463 1F (RQ 6GEeL RIND)? GO TO 50 68.1205 
000246 RA ® KBO) 
000250 GO Tn 448 
c 
000256 So wWRITEILUN,2) KAHO 68.1205 
000256 CALL AMIT(oG9e3m eKARL) 6900709 
00026) KARO 1}) B3rC 6920709 
000263 KARO (4) BROATE 69.0709 
900264 WRITE ¢LUP 92) KARL 6920709 
000272 60 Tn 30 6821205 
€c 
c 
000273 60 EWN FILE LUN & REWIND LUN 
000277 CaLt tal 
Cc 
Cc 
€ 
000300 Bud WHITECGLo1) (MESG(Me2) Kul ed) 
060312 Go To 64 
¢ 
000313 Aye WRITELGLeD) (MESG(Ke3) Ke) 94) 
000325 GO Th Blo 
c 
000326 8U6 Keekac? 
000330 KG0s4 
000331 Gy To 112 
c 
000332 AOS BRITE(GLS1) (MESGIKoG) sKales) 
0003464 < Go To 106 
000345 ALG WHITE COLoLD (MESG(Kel) oKedea) 3 
€ 
000357 Ajg CALL EKIT 
000360 END COMPRKG 
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COUNOUT 
a” 
SUBROUTINE COUN OUT (LINES) on.0209 
e SOURCE DATE 66-0209 COnvERT TO COC 6000 
c SOUMCE OATE 6720601 BRAND NE® CODE 
€ 
c SKIPS PAGE AND COUNTS LINES FOR PROGRAMS BHICH WRITE ouTPuT 
bd € TyPrcat USE IS TO CaLL COUNOUT (WN) BEFORE WRITING N LINES OF OUTeUT 
Pa c Ip N IS NEGATIVE PAGE SKIP 3S SUPPRESSED © THIS MIGHT BE USED WHEN 
¢ You Want NEAT INPUT TO APPEAR ON NEW PAGEs BUT IF TMERE IS 
¢ NO MORE INPUT» YOU wanT TO aVOID THE EXTRA PaGEs 
C 
CCOMPKG,BASCON 
000003 COMMON /BASCON/ 
1 — PHOGRM, KPAGEs LINEs TOF Have RUNG RUNIN (606 
2  SGe FLAGs DOATEs MAKTMs MAXPGe MAKLNG 
4 3. -KLASSe KGROUP, KKUNITS.  KSTINTs  KOORDs IF EOF 
4 & eNOFILE 
; ¢ 
¥ ¢ 
‘ 000003 IF(LINEOCLINES GT eMARLND CALL WEAD(2) 
’ @0001) LINESLINESTAUSILINES? 
; 006013 RETUAN 
4 @e0014¢ ENO COUNOUT 
i ¢ 
z . 
{ 
* 
F 
a 
a e 
Fi 
He ~* 
A: 4-19 
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Spt Ne tome 


000005 
000006 


000006 
000015 
000016 
000024 
000030 
000038 
000037 
000040 
000062 
000052 
000054 
000055 
000087 
0000S? 


SUBROUTINE CROSS (Ae80C) 
SOURCE OATE OF0071) BRAND NEW COOE 


RETURNS Co THE CROSS PRODUCT 4 Kk B 
FOR ENTRY CROS]e RETURNS C UNITIZED 


DIMENSION &(3) 0813) 0C(3? 


«R0$Se) 
Go To 10 


ENTRY CAHOS] 

KROSS#0 

C1 94(2) 9B (3) 0013) B62) 
C2aa( 32061) 0a02) 083) 
C3) eA(1) 0812) 0412) 0801) 
ct2)sc2 

ctlieci . 
TF (KROSS) RETURN 

C1 sSORT (C1 0020020020013) 092) 
Celpeceiysc) 
ct2)0ci2)/C) 
ct3)0ec(3)/C) 

RETURN 

Eno CROSS 


5 iA ae Mev 


ai RENO 


VER Lel 


000003 
000004 
oocees 


000006 
0000206 


OaTer 


eee0oce §=6oe0ele 


04012405060000000001 


7660 
01000 


L) 
S160¢ecee? « 


24113560000100000007 ° 


51580 
10655 


eeee?y ¢ 
0305000005 « 


020000000) © 


034911 


DATEF 


06/30/68 


1OENT oarer 
PROGRAM LENGTH 


SLOCKS 

PROsRAMe® §=LOCAL 

ENTRY POINTS 
ee0ces DATEr 

EXTERNAL SYMQOLS 


cec 
e SOURCE DATE 66.6208 CONVERT 70 COC 6400 
e RETURNS DATE IN 10 CHARACTERS 
enray pater 
vFo 2O/SLOATEF e3072 
oater e8sz 1 
He 80 
sae aSANng 
RJ excPc 
VFO LOZICTIME LAL 91/1047 Q01 271 0O/00 1 B/ANS 
Sas eSANS 
za KSr® 
Oke us 
Pd oater 
Eno 


UNUSED STORAGE lg STATEMENTS 3 SymBOLS 
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DISCON 
@ 
a 
SUBROUTINE CISCON (GUAL oGhOm GLAST snUNDASDTeDTE) 2 
e i SOURCE DATE 6960008 BRAND NEw CULEs REeWI TIEN TO ULD SPECS 
2 c Suurre DATE ©460328 SIMPLIFY GUAL/ULAST/ONU® TEST 
e SOURCE LATE 65.02U9 CONVERT TY COC 6600 
j Cc SOURCE LATE 0600826 
¢ 63 
c ADUUSTS Time STEP TO ITERATE BACK TOWARL » DISCONTINUITY 
c as 
ry c eda - NUMERICAL VALUE OF CRITERION AT eHICH INTEGRATION ID 06 
c TO BE INTERMUPTED “7? 
‘. € Show - PRESENT VALVE OF GOaL 48 
c RLAST © PmEVIJOULS VALUE UF GOaL 09 
Cc RONDA ° CUNMENT COnNUITION OF LIEHATIONs SET BY DISCON TO « 50 
c POSITIVE © ITERATION IN PRUGHESS*® RONDA & STEP NOe 5} 
¢ NEGATIVE = ITkmaTION COMPLETED ImIs TIME STEP $2 
¢c ZERU = NO ITERATION IN PMOGKESS (MUST BE INITIAL 53 
¢ VALUE OF KONOX SUPPLICL TO ULSCUN} S6 
¢ im = LAST acluat TIME STEP $5 
: ¢ OTE ° TIME STEP ESTIMATED TU REACH UISCONTINUITY b6 
‘¥ c 57 
a € 60 
1 Q00011 OGeAHS (GOSL°GNOw) 
000013 TLANCE) o£@S © AMAA) IGUAL ele?) 
000017 TF UKONDA) 294610 
é 000020 2 “ONOxs0 
k 000021 &@ IF (DGeLTeTLANC) GO TO 12 
} 000026 TE C(GOAL -GNGw) © (QUAL GLAST) 261,00) 60 10 8 
2 ¥ 000031 6 KONOXSKOMDASI 
; a 000032 OIE SN TOFU1V 4 (GOAL <GNOm) / (GNGW@GLAST?) 
en 000045 TF (DYE EQ eOe) DIEZIU.POT 
: 000007 8 GLASTSGNOW 
- 000050 RETUGR 
_ ¢c 
a 000051 10 TF (UGeGTeTLRNC) GU TO 6 
000v55 12 KONDxz-] 
000056 Gv Ta A 
q ¢ 
ag 000057 ENN NTScuN 
i 
rt 
4 
si r 
x 
“sg e 
e 


OS aa 


7 DNSITY 


e* 
FUNCTION ONSITY(S) 7 
c SOURCE OATE 68-0209 USE aLTF FOR ALTITUDE 
4 SOURCE DATE 64-0728 
¢ 78 
¢ RETURNS ATMOSPHERIC DENSITY FOR STATE VECTOR $ 7 
e € 60 
€ WRITTEN 7720766 ai 
3 = ¢C a2 
i CCOMPKE »BASCON 
eeeec3 COMMON /BASCON/ 
1 PROGRAM, KPAGE® LINEs TOF DAYs RUNe RUNID (Ole 
2 MSG FLAG, Dates MARTMs AXP, MARLNG 
3 KLaSSe KeRQUP, KUNITS. KSTINT, KOORD> IFEOF 
@ 4 & oNOFILE 
5 ¢€ 
3 CCOMPKG »CONCON 
#4 eecees COMMON /CONCON/ 
xX 1 Pte SRD» SLVs sw, SKPo Reocoy, 
. 2 Gacce 6CON, WBODYs RMOZROs TwOP Is warPt 
\ c 
= DIMENSION 2(5)0A(S) 0815) 
4 Data ((ZUT) e702 9S) 80} 00 0500000100000 540000091000000) e? 
P DATA ((ACL) oTabeSd ah cedeel 08004620601 50678735» o3S487133) ae 
€ 3 DATA ((OCT) 010) 0S) 800%200030160%,00012 7052000156) 1920001185) ee 
* ¢€ 
rd 00003 WeALTF(S) 60.0200 
@ 000008 «6 CLF (KUNTTS.NELO) Wenssne 6 
e00e0i1 on 7 19105 eo? 
00ee13 IF (Mok To2(2)) GO TO @ oe 
? CONT JNUE td 
a ONSITY#0. lee 
z RETURN 01 
| e ONSITV@RHOZROSA (7) EXP (815) OH) 102 
e ae TURN 103 
ENO ONSITY 
; e 
: e 
ri 
= 
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; A 
FUNCTION NOT (Aed) 108 
; c SOUACE OATE 66.0101 GRAND NEW CODE ‘ 
- ro 06 
: @ rd RETURNS THE DOT PRODUCT OF VECTOR A WITH VECTOR B 107 
c 108 
=. id 600004 DIMENSION A(3)08(3) 109 
1 d 
} i = potso. 110 
| M0 3} Iele3 Mk 
: 1: DOTSA LTD ORLS) DOT 112 
RETURN 433 
END oot 
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DTIMEF 


08/30/6n 
TOENT OTIMEF 
PROGRAM LENGTH 
BLOCKS 
paceRame §=LOCAL 
ENTAY POINTS 
000003 OT MEF 
EXTERNAL SYMAOLS 
cec 
® SOURCE DATE 68.0209 CONVERT TO COC 6400 
e RETURNS TIMESOF-DAY IN 10 CHARACTERS 
ENTRY OT IMEF 
VFo JO/OLOTIMEF 924/} 
OTIMES essz 1 
Sxo 80 
sae eSANS 
QJ excPe 
VFO LO/ICTIM01/101/154/9012/296700 10/ANS 
Sas sSans 
zm Li thd 
aXe xs 
Nal OTINEr 
END 
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an 
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Pl 
FyNCTION ELF ta) $34 
¢ source Date 062020) grand NEW cobDe 
: c 13 
Ye c RETURNS THE ELEVATION ANWeLe OF vector & 140 
i ° c 141 
c ANGLE MEASURED FROM AeY PLANE® POSITIVE rowarRo 1eAktS 
e c 11S valve tS BETWEEN oP an Pl. 
¢ 
000003 DIMENSLON 043) 9Bt39 3492 
c 
9e0003 ec2ye0at?) 143 
BcBNeSORT (ALL IOALID OACZIOACEI? 166 
EL Fede 
YF vedas 6£Q2000AND0B 12) 060000) RETURN 
E_peatanz(at2) oot?) 165 
ReTuRN 146 
eno «ELF 
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SUAROUTING EULANG (MONE e ARES ANGLES SNNDERI ¥) 108 
. SHUXTE DATE 620209 CONVERT TO CUC 6400 
c SOURCE DATE h600921 REWALTTEN TO aLlLow ANGULAR RATES 
Cc 109 
= c PHODUCES EULER ANGLES FROM DIRECTION COSINES, OR VICE VERSE 150 
1$1 
e C rOvE - IF POSITIVEs EULER ANGLES AME RETURNED FROM OIRECe 182 
c TION COSINESs [F NEGATIVEs DIRECTION COSINES aRE 183 
c RETURNED FROM EULER ANGLES. 156 
c AAES = ORECTION COSINE aARaAY 155 
c ANGLES © EULER ANGLE VECTOR 156 
¢ NUEKTV = MIGHEST ORDER OF TIME NERI VAaTIVE INCLUDED 15? 
c 154 
Cc 160 
CCOmMPKG,CONCON 
000006 COMMUN /CONCON/ 
} Ple SRUe SLVe SMF, SKPe RB0OVe 
2 Gacle GCUON» eeooy, RHOZRNO, TwOP ls warey 
€ 
000006 OUMENSION EN(10) GAAES (2003) eANGLESI 10) eSNFNI 4) gCSFNIO) 
000006 NIMENSTON SNENO (6) 0CSFND (4) oe SNFNDO (4) e CSFNOD (6) eANGLS (10) 
000006 EQUIVALENCE 


J (SNASSNFN(2)) 6 (SND) SNEND (2) do (SNADM oe SNFNDD (2) ) 

20 (SNH eo SNE N (3) ) © CSNANM > SNFNO (3) ) » (SNBOMe SNFNDO (3)? 

30 (SNG*SNFN(8)) © (SNGD* SNFND (43) 9 (SNGONe SNENOD (4) 3 

wo (CSAeCSEN(2)) ¢ (CSADsCSFND (2) } 6 (CSADDeCSFNOD (2) ) 

@ 5s (CSHeOCSFN(3)) 6 (CSHOeCSFND 13) ) ¢ ICSBDDECSFNDD 1 3)) 

6+ (CSGOCSFN(4)} 6 (CSGDeCSFND (6) ) 6 (CSGONeCSFNDO (4) ) 
000006 EQuI VALENCE 

) (AS @NGLS (2) ) © (AD e ANGLS(5)) » (AON eANGLS (8)) 

2+ (AeA@NGLS(3)) > (BOeANGLS (6) ) 0 (BOD saNGLS (9)? 

3s (Ge A@NGL514)) ¢ (GDe ANGLS17)) 2 {GOD 2 aNGLS()0)) 
000006 EQUIVALENCE 

1 CENKOFNG2) 1 0 CENKOSEN(§)) e (ENADDCEN(B) ) 

2s (ENV OEN(3)) © (ENYNOEN(G) De (LENYADCENIO) ) 

de TENZSENIO)) © (ENZDOENI7)) e (ENZDDDENIS0)) 


€ 
¢ 
000006 {F (MUUVE cL TO! GO TO 30 
an0n07 ONGLES( 3d @ ARES (4 0)) 
000010 CSA ® AMIN CAKES (COeZ ole 
060016 « # ACOS(CS#) 
600017 5NB ® SIN(H) 
60002) 6NABb ULV 19 dE 1303) /SNR) 6200209 
000032 FNYsFULY (ARES (263) /SN8) 6a,0209 
000061 IF (CSrHet Tepe) GO TO 6@ 
000944 FNK @ ALES C201) 
ey e 000005 PNY w AKES(3e1) 
ey 000ne? ao NZ © Oe 
% 000050 Siva B ENR 
t ° 000082 CSA 8 ENY 
e 000053 a ® ATANZ(COAe DNA) 
i} 000056 G @ SCPALENCZ) oAMES(Z0TIF 
&. 000ns6 TF CAKES (Ooh) oh Te.) GOeG 
4 000071 Ite 4 
; e 000072 TEINPERT Vet Tod) GO TO 28 
¢ 
9 o00nTs HOSFULV COMXES (70 31 /SNH) 66.0209 
000103 TE USN eEWeOe) BL BSIGNISQNT(AXES( 701) OO 2AKES (722) 002) »AKES (OS) 3) 
ooo117 SNAD ® CSHOHD 
600121 ENXD © FOTV Ie (AXES (693) @ENASSNBD) /SNA) 62.0209 
000132 ENYD # POLTVE (AXES (553) ENA@SNGO) /SNA) 68.0200 
000143 €NZD ® O6 
000106 an © FOLVGENYOSCSa) 68.0209 
000181 TF(CSAsEQe00) ADMOF OT VIENKO/SNAL 6800200 
00016) CSGD & DOTLEN(2) oAKES(Se2))¢ DOTLEN(S) eARES(202)) 
000201 GO ® FUTV(eCSGD/SNG) 68.0209 
900205 TF USNGsEQ.0e) GDBSORTIENXDE*2cENYDE*O2) - : oa = 
000215 rs? 
000716 TFINMERTV.LTe2) GO TO 28 
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ie ¢ 
£ 000223 400 w FOL Wt 1aXtS(1003)-8O®SNBD) /SNB) 68.0208 
Hy 000236 SNBUD & <SNKORD®#2¢CS8°H00 
4 0e0zel ENRDO ® FT ¥ 1 = (ARES1993) 92s CENADOENBD~ENKOSNAOD) 75NB) 60209 
‘ 000256 ENYOD & Firlv(o (ARES ( R93) 2. CENYDOSNBN~ENKESNBOD) /SNB) 68.0209 
080270 ENZDO * 02 
000271 and « FOLVILENYDDsSNs®aDee2) /CSA) 68-0209 
000300 GOD w FULV © (DOTLEN(2) ARES (B02) ) e2e@00T(EN(5) eAXES (502) ) 60.0209 
x * DOTLEN(A) eARES (202) ) =GNOSNGD) /SNG! 6000209 
000337 tr # to 
e? C 
000300 28 DO 20 T82e11 
000302 29 ONGLES(I? @ AaNGLSE]) 
090307 WE TURN 
c 
Cc 
000350 30 OXESt]e1) & ANGLESI1) 
000354 DA 32 Jeers 
° 000355 SNFN(J) & SIN(ANGLES(J)) 
000364 32 CSFN(J) © COSC(ANGLES (JI) 
* 900376 IFINDERIVeLTe1) GO TO 60 
090000 00 Je 58206 
890402 SNFND (J) © CSFEN( U1 ANGLES (293) 
600006 34 CSFNO( J) & @SNFN( J) @ANGLES (03) 
Qooars TR INDERIV CT e2) GO 10 40 
Go0ate NO 36 J8Zee 
000020 SNFNDD (2) & -SNFNIS) CANGLES ( Jo 31 ©@29CSFNI I) OANGLES (Jeo) O02 
000031 36 CSFNAO CIE = @CSFEN( UZ) ANGLES ( Jo 3) C8 2@SNFNI SL) OBNGLES ICG) COA 
€ 
On0ee4 20 AMES(2e1}e CSACCSGHSNAOCSHESNG 179 
90005¢ AXES (301) = SNAPCSGOCSASCSBOSNG 180 
000056 ARES(Ee)) a SNRESNG ii 
000057 ARES (292) B=CSAP®SNG~SNASCSBOCSG 182 
boose2 ARES (342) BeSNAOSNGOCSASCSBOCSG 183 
000465 axESlae2)@ SNA®CSG yaa 
0900066 AKES(2¢3) 6 SNASSNAR 185 
000470 AXES (393) S=COAPSNA 184 
QoGaTL AaFSlee3)3CSh 187 
000073 TE CNUERE VL Tel) METURN 
c 
@ 000076 Pe = CSADOSNG*CSROSNGD 
dO0sS02 VY = CSADECSU*CSROCSGD 
000605 OXESISel) 3 CSAN®CSGeCSA®CSGD-SNADECSAPSNGHSNAOPR 
000515 AXES (pel? @ SMAD®CSG¢SNASCSGD-CSANOCSHOESNG@CSAOPK 
00n825 ARES( Tel) & SNHD®SNGeSNH®SNGD 
0003830 AKES(S02) & —CSADOSNG=CSAPSNGN=SNAD*®CSB°CSG-SNASPY 
00054} AXES (602) © =SNADSSNG@SNAPSNGD= CSADOCSBOCSGeCSACPY 
000550 AKES(Te2) = SNUD®CSGe SNBOCSGD 
9005546 AXES (5¢e3) © SNAD®SNBe SNASSNGO 
000557 AXES (603) @ =CSAD@SNB@CSASNAD 
000862 AX€S(793? © CSbD 
000503 TE(NDERTV.LT.2) RETURN 
¢ 
0005866 PKO = CSBDD@SNG « 2,eCSB8DeSNGN +» CSBeSNGOD 
000575 PYD & CSAND®CSG * 2,.eCSBO*CSGN ¢ CSBeCSEDD 
000603 AXES(OBeL) @ CSANDPCSGe2.*CSan*CSGDecSaCS$G00 
ee x @SNAN0 OC SBOSNG<2. SSNADOPX=SNAEPXO 
ie 000622 AXES(090)) = SNAND®CSG22,#SNAD®CSGD*SNAPC§GDD 
: ‘ *CSADUPCSBOSNGS2. PCSADEPKOCSACPAO 
“a 00060) QXES(IO01) & SNBODESNGe2, @SNBD®SNGD + SNBOSNGOD 
ot 0004650 AXES(0B02) & =CSADDOSNG@2 0 ®CSAD®SNGD-CSASSNGOD 
° . x *SNADD®CSRECSG 2c SSNADOPY @SNASPYH 
: 000667 {thee hese & -SNAQDSNG=2 o @SNAD®SNGD~SNASSNGDD 
i ° *CSADLDECSB®CSG42 0 *CSADEPY eCSaewyD 
34 000706 ” huest yond & SNBOD®CSG02. @SNBNECSGD*SNBOC S5G00 
000735 ARES(OBl3) & SNAQDPSNGS2.*SNADSSNRD® SNAPSNBDD 
800723 AXES (0903) © eCSANDOSNB~2.*CSAD9SNBO~CSAC SNBND 
000732 AXES(1003) © CS800 
‘é Ooor36 RETURN 
Gq0736 END EULONG 
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cE Toent Folv 
: - 000005 PROGRAM LENGTH 
4 e @LUCKS 
4 : voouno 000005 PROGRAM® LocaL 


“ ENTRY POINTS 
000001 FOIV 


SOURCE DATE 68.0209 CONVERT TO Cnc 6000 
& FUNCTION CALLED WITH ONE aRGUMENT = FOTY RETURNS Twat aRG 
UNCHONGEDs UNLESS IT IS ILLEGaL (IE = INFINITE OR INDEF INT TED 
In WHICH CASE IT RETURNS ZERO, TYPICAL USE IS eee 
INSTEAD OF keA/B 
WAITE ReFOiviasad) IF B MIGHT BE ZERO 


7 de eames - 
e@es0oe8 0 


; ENTRY FOIV 

: 000000 dans1i2ooa0e000ccen? VFO 2e/OCFNIVs 36/1 

ee 00000) Folv assZ 1 

e 000002 Sain Sai al 

10611 8X6 x1 

; @ 0351000006 ¢ oR K1eSETZRO 
000003 0361000001 ¢ oF ADOFOIV 
000006 76600 SETZARO $x6 80 

0200000001 « wp FOolV 
! 900005 END 
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SURMOLTINE FLAC 


SPurce 
Source 
Source 


Source 
Source 
SuUNCE 
SuuRcE 
SOURCE 
SuuRCE 
Suurce 
SOURCE 


SGURCE 


vale 69206013 
UATE 6920325 
UATE €0e}265 


VATE o@e0T71T 
VATE 66-0311 
Date 68.0209 
DATE 68-0208 
DATE ©8e02u7 
DaTt 6760)1 
DATE 67.0726 
OATE 6700707 
67-0707 
Lal— eFefel2 


(SeCeCr) 


CORRECT CALC Ub aL Prac 
WEeARKANGE SE UP OF ANCe DUMMY (e"6) 2 AND OD ARRAY 
USE KUNVERG INSTEAD UF LIVERGE 


ALSU MAKE AMBIENT OPHESSUNES TMRUST©CURKECTION 


UEPEND UN SUNIC SPEED aS WELL AS OENSITY 
uSe QUaURATIC FOWM FOR ALF UF CNL eCN2 
PASS CP (aLSO) TU CAaxX1aiL AND CNORMEL 


CONVERT TU CUC 6400 


INCLUDE DRAG BREN FINDING ALPRA 
CORRECT ALF FUR NO@MNVK, NOPRESP FLIGHT 


AMFND FORMULA FOR ACCEL wI TH RESPONSE 
CALL TRAIOERK IF ERROR CONDITION 

CALL CaATALe CNOKML FOR AERO COEFFSs ALSO 
ANC IS IN RaviaNs IF *FLAGe IS NON@BLANK 


ADAPTED FRUM FLACUN 


CUMPLTES aCCELE RATION OF HOCKET #1TH UVAMPEL SECOND-URDER 
Oe SPGNSE TO NORMAL ACCELERATION CUMMANULS 
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ce 


- CURKECT ROCKET STATE YECTOH, FULLOWED BY CURRENT 
ROCKET ANGLE @OFmaTTACK VECTOR 
© ChARACTERISTICS VECTUR FUR MOCKET TYPE 


C(RSTAGE ee}? 
C (RMSTAGE 92) 
C(RSTAGE» 3) 
C(RSTAGE 94) 
CURSTAGE 5) 
CiRSTAGE 06) 
CIRSTAGE 97) 
CURSTAGE 08) 
C (KS TAGE 09) 
CRS TAGES 10) 
CURSTAGE 011) 
C(RSTAGE 012) 
C(KSTAGE 613) 
CURS TAGE 916) 
CUKSTAGE 0} 5) 
C(KSTAGE 016) 


VACUUM THRUST AT IGNITION 

NOZZLE EXI? AWEA 

TOTAL MASS AT STAGE IGNITION 

TOTAL MASS aT STAGE BURNOUT 

STAGE BURNING TIME 

MEFERENCE AKEA 

SUBSONIC aAlTal FURCE COEF, 

SUPERSONIC AAITAL FORCE PMT a 
SUPERSONIC AXIAL FORCE PMTR B 
MYPERSONIC AATAL FORCE COEF. 

MAXIMUM ALLUMABLE NOKMAL ACCELERATION 
MAATMUM ALLUMABLE ANGLE OF ATTACK 
NORMAL FORCE PARAMETER CN) 

NOKMAL FORCE PARAMETER CN2 

KESFONSE DAMPING FACTOR ZETA 

RESPONSE UNULAMPED RESONANCE FREQUENCY 
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CURKHENT STAGE NUMBER KSTAGE 
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RPHASEZ©] — PRE-EUKNING COAST 
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STAGE 1 IGNITION TIMe 
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STAGE 3 ToONTTION Tee 
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C 
c NeGAlIWE © EXPOFENTIAL BURNING KASE* CONSTANT 
c THrUST©TOeMASS HATIO 229 
c 230 
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c We TURAEL IN S16) °S8(20) FOR INTEGRATION FY STINT ALONG WITH Se 2eu 
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J VUC9) G13) 0003) ewORKK (mR) 68.3205 
RECT VALERCE Cid yd oP Ty de (902) CAE) oe (D693) OAM) 9 Mile) Amare (UISD oTede 
V teed eSmeR oO (DIF) CAL) (Ue) oA OLLI) obra (LITC) CA) 60.066) 
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Testy) 256 
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Ir (nSTERePE eu) GLU TO }1 eos 
Viele .Seit 256 
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Cunt ve 260 
Ke13 2e) 
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000 Te meu mw PASTTY OS) 66-060) 
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900165 ls Files, ee? 
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000171 CALL GRAV(SS8G? 290 
000176 VRRMAG(S{(S)) 
2 0002u3 er 
t 000ALS wSOmat@F) PV ISHEFZAMASS) 68.9209 
000213 GSA MaSO" . s 
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Qunesa2 CaLe VECLIN CARDS (5) pLeeGed(6))} 3 
000ceC me TUM. H 
a Cc 4 
ONH2HI 2.4 BSSTenN 215 TO RUM paAK 68.n20b 3 
O0N2bs BLE soL Pram 4 
7 000244 GL TA MARECA ‘ 
000273 415 ASSTat 21e TO KUeB aR 64 eM20h 
=: 000276 Ge TA MAKECA 6B.0c0R 
00030nN 235 GLIMAWOSCSUBROUSAM 
000352 GEIMTPRBGL IF AROO@CUS (ALF) © (F UmM@QSUMOCA) OSINIALED 68.020R 


000316 Te (ANCMAG CGT 20.) GO TO 210 


a et a 


nud3e? 


00U3¢3 
00N 324 


O00 346 
06035? 
000333 
O0U334 
000337 
QO 300 
000347 
000364 
000351 
Que 456 
OO ie 
00U 3e7 
yooa7t 


oCyus7? 
O0ON 376 
OUGa:? 
OOOaits 
O0He0) 
N00eV3 
0004S 
ddGect 
O0Q817 
000420 
0006¢3 
0008 36 


O0hees 
00y4>5 2 
OL0as4s 
OI04e% 
00na7s 
QOuaTA 
VdUsrs 
O0NKue 


GOOseT 
ounsid 
OUugSlA 
000517 
acgSes 
000539 


070532 
000536 
000535 
00056? 
000545 
0005446 
000992 
000554 
0u055? 
0005e9 
000562 
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Fue 


gly 


ee 


2% 


alLeaan, 


aSSqet 233 | AUmOAR 
5U To wAKus 


LCE Aah} 

bub Bauer be 

1h (Gt TMArL eb Taahh) GO TO €20 

CGUAL BANM/ USAR 

BoSian 216 TO numban 

ORFRRELONS 

CALL RONVSET (oC) eo URE e ORK) 
ALFREL IV ICOCAL/UCUA) 

bu TO wanech 

Tests SurneCGual 

Lb CRONVENG CUTEST OAL E omORK)) 213993990620 
CALL TRA DERRI2@raLHnMk TROUBLE IN PLAC @) 
HE TUR. 


ALPmMwasALF 
IF (hUCAGetweOL NK? GU TO 222 
ALF BARCMAG 
ou Ta 238 
AL’ eaALPHary 
16 (Gt LMT ee eb TeANCMAG) GU TO 230 
ASS1on 226 1G RUMBAR 
CIMYSEC Tv (aNCmaGs LOSNM*FOTVIFOM/UCLm))) 
On€ wma, 00} 
CALL RONVSET (ou) eURE HH OM OMK) 
ALE sete T¥ICTMYsUCUA) 
ThiCt ee Mb eVe) SLE # 
LeCN] eSERTICN1 @O2 eb e@CNZOCTRY) I /(20°CN2) 
GN ole MARECK 
Asstt. 22e 10 cUMBAR 
SU TO wAKECA 
TEST EC SUHN®USNM®CUS (ALF) © (FOM@USOMECAT MSIN CALE IM ANCMAG 
ASSTGn 224 Tu RLPBAK 
TF (RONVERG (TEST oALF ewORK)) 22593990236 
CALL Tha LERHI2emMaALRPRaC TROUBLE Ip PLAC @) 
Re Tur 


ALRMACEALF 

ALF AaaMIRy (ALPRAC PALPHMA s&LEHAM) 
ASSTG:. 2322 TU KUMGAR 

GU To BARBUS 

CSABCUStaLb rad 

SNHESINGALERA) 


ALF SALPRA 

ASSIA 2ie TC KUMMAR 
GG TO wanercCa 
BASF ir eCAeUSCi 
ASSInt. 230 TO RUMBAR 
GO Tr wARECK 
ANECSLANO@USNE 

IF (ANCMAR eGT ele) GO TY 26b 
CUMLsUMEG LEP? 

CUN2Z a? oe PLUME GACLETA 
Ou 24 Naie3 


66.1265 
66.0204 


66.1265 
6620209 


O7eNTHT 
HTN? 
OTeNTL? 
67.0707 


bkenheny 
04.06) 3 
burleus 
on neny 
67,0710 
67.0717 


O5eN204K 
CHeNetH 
66eNEUA 
66.020H 
6b.1end 
Oh.0eL9 


Fe 


Ep tos 


00h? Simer 77 SCULL@ EVE (ND Ut) POCUN2#S (NO le) 

aves??? Otte pabb be if IND SURLG) ra Pair an) 
OOVS7T 23m Str e7) BAAPLIA) Gila) 

O0UELO rr Tut. 


000410 eeu Ir (OskR Gace ede} GU Th 36 

60061) LU Je hBLe 322 

000A) 4 cry SIN 7) BAKO (WL Civ) OCSACACU IN) SNA) CANE LACU IN) PCSARVU(N) OSNA) GIN?) 323 

0006 32 Re TURN 326 
id 

000632 je CSF BCCSLALFA) 325 

000834 Sng zSIn (arf ad 326 

000640 Cun) sUrEGaert? 327 

ondsel CON2 8? ePGMEGASZETA 32a 

Vo0ee3 CUNZct OI viansSna) 6.0269 

000656 tO ja NB] O93 330 

000653 Steer 7) SCUMLO ICSE OU LN) O StF PACUIN) UIN)) CONZ®S (No) ®) 33) 

00nhHA Cth ef OL v (O(N) JUMAG) 64.0204 

000673 3” S (he 7) BAA®U(N) ©CON3ZO1O(N) @CSACVU (ND) OGIN) 332 

yoo71) ME TURT 333 
¢ 

0007211 399 GO J MAKECK bbe eb5 
C MARECIN 

000714 ey CSUnvaC iC el CvmMaCMe ALF oC eKSTAGE e CP) 64.0311 

o007eas GU TO RuUrdaR 67.0707 
¢ 
CN ARECA 

0000/3" OP CABCAATAL (VM ACR OALE 0 C oKSTAGE 0CP) 68.031) 

0007TH#4 9G Tr KRUPRA~AK 67.0707 
C ‘eRUYS 

UdO?a? 446u UU eae NELOD 

060764 VUCN) SFOLV tS tivea)/v) 68.0204 

600755 O(N eS (help) evUtr) 316 

O0U?A) eee ACUI) ab TV (ANC (IR) JANCMAG) - 66.0209 

900773 Gh aosamaG (ud 67.0831 

OCA? ?a ALPMASALFA 

000777 IF (UMEGM ove e Oe) ALPRARSEPA(USVU) 

0010045 GU TO KUMnAR 


oo}7i3 End beL aC 


FLIER 
SURMOUTING FLIER (Se CR eKONDASDTE} 335 
C SuurRce VATE 6601205 EXTEND REVISIUN GF 68.050? 
Cc Suomch Pole nBetb]2 OUNe1 JUMP iu ST, 6 Ih ANC ubper INED 
Cc Scurct balk 6420507 END MANEUVER STAGES ON ZERU IN Ch aRRay 
Pas fa Suveee UsTe 672023 FIX FOR ROUNL EARTH 
ig Suumre DATE AhelOv) FAPANDED FOAL CAFABILITY ; 
Cc 336 i 
Cc CeNtraTeS MAKLUVEW COMMANDS FOR FLIGHT Balm rr Oormam 337 
C 33 
c S = vericre STATE veCTOn 339 { 
C ce = FARTICULAN VEHICLE CrAHACTEWISTICS 340 
= € ce lpd=chty2) - SEE #FL ace FOM DEFINITIONS 34) 
2 Cr V1L3}9CR(15) |= MANEUVER STATUS INLICATORS 342 
a ¢ ce¢13) @ NUMBER OF CURRENT MANEUVER (ZERO TE 343 
c KU FLIGHT bMOGHAM 15 IN USE) 34% 
¢ Chile) = VALUE UF MANEUVER GUAL AS LAST TIME 345 
¢ STtP KEGAN 346 
c CP OTS) - TIME al wrHICe CURRENT MANEUVER BEGAN 347 
C Crh ledbeCR (20) = SPECIFICATION FOR FIRST MANEUVER 3a 
¢ In FLIonl PRUGRAM 349 
c cr ()A) | MANEUVER ACCELERATION MAGNI TUUE 356 
f ig Creln) ~ MANEUVER DIRECTION (MEASURED FROM 3o1 
i ¢ HORTZONTaLs AKULT VELUCITY VECTOR) Be 
Ey c ZERO IS mom LZUNTAL TO PRE HIORTs 
m c 9 DEG 19 UUme 
ae ¢ ch(LR) © FATE OF CHANGE UF PANEUVER CLRECTION 343 
¢ Cktiy9) - MANEUVER GUAL (VALLE OF CKIFEION aT 356 
: t waiCh maAncuvER 1S 10 END) 355 
g ¢ cP (20) = COVE DEF imING GuAL (SEE FOAL FOR UEFe) 3%6 
: c Cel2tecr (25) © MANEUVER 2 LEFINITION 35? 
: c tlCe FICe 356 
L z fa itt - €STImaTen Time STEP TU BRING VERICLE TU GUAL 389 
; d C FUNDA © CURRENT TTERATION STAIUS IN ENDING MANEUVER 30 . 
3 (SEE #HTSCON#® FOR DEFINITIONS GF OTE AND KONDA) 361 
c 3e2 
CC MP RG eg hASCON 
OGOnnD? Creme Sen SC Us 
1 PIVUIGHE RPALE » LINE. TithUAYe RUN, MUNITD (0) 6 
2 MOGs FLAG? Dal&e MaATMe MAAPGs MAALNO 
i 3 KLASDe KGMOUP » KUNTTS¢ KST Tile KUUFDe IF EOF 
‘, @ oNOF ILE 
§ COON RKG aR ASCON 
4 ¢ 09715760 
COO! HROg [TC 
4 000Cu? COMMON STRTCONZA 
, 1 Te RSIEKs NUAFeUe SQANDA (S56) 
j OCH? GIMES STUN DOT eeede = XIMT 1294) : 
: 00UNnT Be VALenCle (S007) eAINTO1) eSOANDK(7)) ‘ 
CC Mw RGs TNTCON i 
Cc 09715749 r 
¢ joe 
a 0000u7 LIMES STUN SUPHD OCH L]1) OM 3) U09) eANC (3) 0713) 67.083 a 
anong7 OutT&A (Zehecets? 07.983 
e voonn7 Etulvabes Ce (Abbe LEU) 66.%60) 
c 
g N000UT wryrMCR (CL 3PPOSeel ic jes 


? ; aduot? 16 (eS tr eftet? GU TO 2 any 


O00) 39 BELSCR (MN ROG) 66.NbuUL 
enor ys Teele ste) GE Tu 12 o8.1eUs 
000014 GUBLECKleNvRe 3) 37u 
OvOu2D IF (LF web GeHrTIME INT) GOLLEGUAL SCP II5) 
vooves CALL UISCUN IGUAL oF OAL (SeChk (MNV ROG) OCH Clad enti KeL TOUTE? 
O00) 1F (KONDAcGE CU) GO TO 2 373 
0000465 W NVA SP NY OS Ble 
0000467 CPCL UZ tMNVROE LL /S 
000053 TF (CR UMNVHOO) oEleebe) GU TU 1) 66.12U5 
000085 Ch (La) BFODL (Se Ch (MNVR Oe) ) 
00053 CPhcisissiyd 377 
000%+S P TE ACH UME VR) EUs Ce) GO TO LE belles 
0000H? Le (ROORUcEe2) GO TO 5 67 NHI 
000071 CALL ChUSLIS(5) olen? 67.9623 
O00074& TF CnC) deb uede eANDe WIZ) Eee eANLo MIS) ckWete) milial. 67.9823 
000105 GO Is 7 o7.n8¢3 
000119 CALL CRUST (S(S) Sl 2) en) 6760623 
oou}13a TF CHEL) ot eGe eURe MI2) otibeDe aURe Mi) eMEoUe? GU TU 7 67.023 
O00len CALL CRUS) (20512) om) 67 etleed 
000) 3) wv 13 67.0823 
000) 34 CaLlL CROS) (515) oho) 67,0823 
000)37 POTECh (MNVRO 1) 0151) ) CH] 5)) ®CPRCMNVRO2) 3h8 
00067 CALL VECLIN (CP IMNVRI@COS (PST) oe 601205 
CP IMNVR) @SINIPST) eUe ARC? 6601205 
OOO171 Azle 343 
0co)73 Th CROT( ANC OCH) oh TeQe eANUe KUNDAGtoy) KEM), 68.2205 
000212 OU 9 [#103 346 
o0ncls Cri) sxMancly) 307 
000221 RETUOP 349A 


9n0222 CPL) B06 o8.05C7 
00223 RONDA 80 6621205 
000224 CaLh AMIT(eJeUesCP) oBenole 
000227 RE TU, ov.nele2 
000239 Eon FELTEm 


SURMOUTIAE FLIGHT CAREL eQUIDE sSTeCeCPe ISTEP SGUALFENe VAL sKUDe INSe NST? 
SOURCE DATE 6920909 RETURN UPON KEACHING [NPUT LIM{T NST 

SUURCE DATE 6920709 USF WRITE STATEMENT INSTEAD OF CALL TO wnITIT 
SOURCE LATE 6920328 FIX SAVING OF ACCELERATION VECTOR 

SUURCE DATE 6920109 GRAND NEW COLE 


TWTEGKATES & GUINED TRAJECTORY tu a SPECIFIED ENU CONUITION 


AAEL = SURROUTINE FOR CALCULATING TRAJECTORY ACCELERATIONS 
GUIDE 7 SURMKOUTINE FOR CALCULATING GUIDANCE COMMANDS 
‘ - ARKAY IN WHICH INITIAL TRAJECTORY STATE MUST BE 
FInST ENTHY AND FINAL STATE PRODUCED BY FLIGHT wILL 
wt LAST ENTRY 
ARGUMENTS AF SUBROUTINES AMEL amy GUIDE 
TIME STEP FORM INTEGRAT TUN 
FUNCTION RETURNING VALUE UF CRITERIUN FORK ENDING 
TRAJECTORY INTEGRATION 
RUMERICAL VALUE OF CRITEMION aT wmICH INTEGRATION 
1S TO me TERMINATED 
AKGUMENT FOR GOALFN 
CONTROLS INSERTION OF INTERMEDIATE VECTORS INTO 
ST aRRAY eee 
INS20 GeTS FIRST anu LAST VECTORS UNLY 
In. Saen GETS VECTORS EVERY N Tyme STEPS 
Ingen GETS VECTURS EVERY N TIME STEPS PLUS aT 
FVERY GUIDANCE ANDO STAGING UISCONTINUTTY 
= NUrBER OF FINAL STATE RETURNED IN ameRaY STATES 
ThhCwlant © BEFURE CALLING FLIGnTe SET NST TOU 690099 
THE MIGHEST NUMBER OF VECTORS THE ST 
ARRAY CAN MOLUe FLIGHT wILL PRINT A 
MESSAGE aND HETURN IF YOU HIT THIS LIMIT 6900909 


ac2 
CuPPRKOe TRTCON 


ocnitA Cument, sTt. CONS 
1 Die KSTEP, NOTbEGs SOMNUA (Sa) 
O00 A CIMERSTON SOC}eeG) 9 KINTOL2 06) 
0000146 FwOTVALe*Ce (SOUT) eATNT IL) eSOANDA(7)) 
CCOMMPRGAINICON 
09715769 


¢ 
CCl KG ee MUL CON 
100014 COmENP, seul COW 
1 Seer st) Cee) 
CCCHMPKGeMULCCN 
C 


09715769 
gout LA CIMENSTOR Silte@) 

OOdulLs EuuTvealLErce (SeSPULS) $ UsATA (SReUsy)? 

0000]4 ANIME STUON ST 41 0908ST) eC (3040) 9CPR 120) oSave (396) 


ed INITIALIZE 

O00uU14 CALL AMIT(39DTeSave (393)> 

OlUN22 Dieure se ISTEe 

000030 KoTEUsO 

000031 NOIF EC a] 

000032 TEICELole) NE eUe eOKs Cl201O) eNE Oe othe Cl3elO) eNEoUe) NUTFEQR2 

000066 GUALaVAL & [FCRKUDCEQSBMTIME INT) GUALBGOAL 6ST (101) 

00005? TNastaHs(yhS) 

00UNSs MANS STOT AL IFEW HH eo NGG 


aoocs7 

OO0Un2 

gOOuTS 

4 O0O1os 
“ o0olua 
000107 

000110 

000333 

000117 

000)21 

< 000127 
0u0131 

a 000133 
000] 3e 

OO0141 

000167 

O001646 

000177 


0002v0 
d002v02 
000203 
0002056 
ooceli 
00216 
vooz20 
OOH226 
000241 


000254 
> 000255 


900257 
000276 


MRAM TROT CRORES 


000326 
0003 3+ 


000356 
009357 
000360 
000345) 
q 0003he 
q 000365 
Q00777 
Q006u) 
000405 


s O00eum 
000410 

a 000444 
ON041A 

HHte25 

y0ns 3? 
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2e 
24 


100 


ee 


yue 


ne 


RAL SKAPNSRKASEC 69.070 
GLASZCALe eh (STsKOD) 64.0709 
CALL VISCO (GOAL eOLAS sGLASORXE eDTeL le) 69.0705 
NOTS) 69.0704 
It (RAE eNE CG) GU TO 306 69.0709 
NSTas 69.0709 


OLSTalbeSeStiis)) 

COLL ARTI CL OSNCIFE OST sS) 

LSTesel & (26 

Ib (TSTEPLGTeOe) GU 10 22 

LSTPa-1 > Lays 

ASSIGcR 26 IC Leak $ GU 10 450 

A9Slob Lig Tu Loan 

IF (OT/TLEPTeGTele)d) DTaTLEPT 69,0328 
ChCELS SCHL ZIECHI 31 E06 

Tf (CROUL3) Nb ee) CALL GUILE (SeCPemamel th) 
Cath wkEL (SeCoCP) 

ASSIGN GG TC RUMBAK ¢ GU TO 400 


-—e-—e es se wwe we ewe we wee Fee LivTEGHATE TRAJECTUNY 


i VIMaNTSsrsTer 


Cail STIat 

LAMB Ab 6900324 
IF (KSTEP ott 00) GO TO 302 

SAVE (lel) sS(Hel) & SAVE (2k eSIGe1) S$ SAveEtJol) #St}00)) 

IF (MOIFEGeEGel) GO TO 102 

Save tle2) 291602) S$ SAVECZs2188(Ge2) % SAVE (Ze2) 8511082) 

TE (CP LY 3) eNE Ce) CALL GUILE (SeCPekxMeDTM) 

CALL AXEL (SeC oCP) 

Ie (RSTEP WE SO) GC TO 100 

TLEFTOTLEFT@CT 69.0374 


ed CheCK FO END CONDITIONS 
CaulL CISCUNIGUAL sGOALFN(IS 9KOD) eGLaSeKat sNTeDTE} 
TEL LE L2eANOol GE e7} 

CALL CF TSCON(CHIL) Shel) sGLSTeKASeOTeDTS) 
NSTEPBIS(CLaLIeTFRSISTSTEP © 200} 
TEINSTEP SGT oT 0000 eORe KXESKAMOKAGCOT630) GO TO 500 


IF (RKE) JOUeLUes200 


(F (RAM eivb eo eUWe KKSeNECO? GU TO 19uU 6900324 
DVatster 

IF CUT/STLER Tel Tele? GO TO 206 69603¢H 
Giatiert 69403¢8 
CALL AKEL(SeCoCP) 6960328 


IF (IP Sekued) GO TE 90 
TF (MOU CNSTEPS INA) NEL} GU TO 90 
ASSIFD Yo TO KUMbAK § GO TO #00 


oS 2 8 ee 8 ee 88, eee ee END UF STAGE Ok MANEUVER 
TF CINSeL ECO)? GO TO 2)0 

Swell SSavE (Lee) S$ SiGebdwSAVEC2 0%) S$ Slusel) SSAVE (304) 
ASSICA tle TO KuUMBAK $ GO TO 600 

UIsO 01 @TSTEP 

CALL AXEL (SoCo Ch) 

ASSICh 210 TO KUMBAK $ GO TO 4600 


-=ewr rene nrereenaen te ee FLEw PAST DISCONTINUITY 
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been eee. oe eT OEE ere 


000633 
VONae3 
000467 
anges] 
Q006s2 
000457 
OU04H1 
0006A4% 
OU046h7 
OULa/s6 
ONO47s 
BOHS&H?A 


Onns06 
000533 
000675 


o0nses 
Vnnse7 
von 37 
000545 
000452 
00554 
06045e) 


000462 
00ers 
000574 
000610 
000412 
NOOHA 
000ARer 


ONtimes 
O00A76 
N00ee7 
O0063¢ 
OOnneA 
0004) 
O00K44 


O0nASO 
OOVase 


O00esA 


d000n6n 
000862 


00776 


0007/m 
0007.34 
000736 


oO 


fay 


Ayn 


190 TF CKXS.GT.) eUKe LAMeG1eN) GO TO 200 

SAVE (104) BSAVE C]) S SAVE 12e6) #Save (2) $ SAVE (306) BSavet (3) 
Puu TF(KXS.LT.C) GU TU 4Sn 

TF ckxe.tT.at 68 TU LIA 

StBe 1) sSave lye}) %® S(9el) SSAVE (291) S$ StyuelI BSAVE (30)) 

TF CNP TREC et bet) GU TO 202 

Ste?) BSAVE (102) % SiDe2@SAVE (Zed) $ S(1U92) BSAVE (302) 
PU? TR LTISTEP SLT ¢O.) GU TD 206 

LT mAM{ NL (CTE eUTM es DTS) 

GL TN GH 
av@ DI ZAaAKLIDTEsuTMeuTS) 

Geo Tr 90 


ce tee te ee eee eee eee GET BACK IN PHASE WITh ISTEP 
210 NleleRSe(nSTEPoL @TSTFP = Slie}) 

CALL AXEL (SeCoCP) 

Gu Ir 90 


ee tnt UF TRAGECTORY 
AOSTA Jue TC KLPsAR $ GOO TO 490 

IF (PLE EGeL Ted «OKs [NSeEGe0) GO TU 306 

CALL AMIT(1L0eS (Led) eStle))) 

ASSIcr 3ne TO aUMMAR 89) «GO TO 400 

an@ NSIEPSTH] 

Zuo CALL AMLT(3eSave (103) ebT) 

KE TURIN 


wo 
ze 
Nc 6 


Cl a i el ALNSERT VECTOR INTO ST akRAY 
wud NST= ST] 

CaLh ARIT()OedSeST iL enST)? 

Le UNS TeWE ohaK) GC TO KUMRAK 

CALL COURUUF (3) 

weETs (@eeC1) wT 
SUL FORMAT (SC Jud = FLIGHT HAS COMPUILL@EITE® STATES -#79% ©) 

ie To 302 


-— 2 ef we we wm ewe ere ere eh hl Ohl Fl we le St awCn Um STAGING SIME 
ol LaLeLsTP 

Trrb TSleee 

Th (Leb Te? otRe LeGTeole) GU TU 452 

TLEP TeCRILI<$1)) 

TE CTL EFTSISTEP aL Te eff'l) GU TO 4509 
“52 KASEN 

GO T) Liban 


Ee ee ee ee ed OU A ee ERROM MESSE 
Seo CALL wERI.1) 0) 
CALL STULL (Om FLIGHT RAS FAILED Tu ReaCr a SOLUTION “6 
. LOLI DoxeSeNOIFEQ) 
CaLlL LSKTP (3) 
Cote COURUUT (3) 
AMI TE (605U)) KAMeuTMeCRIL3)» 
. RASeUTSeCHIL) eGLSTe 
° RXt ot TE SGUAI #GLéd 
Gy) FUWMAT(® END MANLY FLAG & #160 vt & @Fe.3e MANEUVER @ SF 6/ 


° ® tw STaGe FLAG & @lae Ut & @Fe.3e6 GOAL @ SF leas 
° e LASTOOAL = 9F11.3/ 

. © FL FLIGHT FLAG & ST4@ ul = ¢F6,3¢ GOAL @ MFileds 
° ° LAaSTOUAL » %F1).3) 

CALL THAI EHH (CIR } 

GO Tr 3na 


Eng) PELIGeT 


69.0328 


69.0328 
69.0328 
6420328 


69.0376 
69.0328 
6920324 
69-0328 
6900326 


6909709 
69.0709 
O90 O7NY 
6900709 
69-0969 


69e032R 


69.032 
69.N3CR 
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ANNAAANANANAANRAAANMNAN NADA ANANDA DNDAN DA ANKRANnRAADDVAAMIANYN BNANnNVANANACNA 


SURMOLTINGE FLINITITLE MONE sCoCP) 
SOURCE DATE 69.0718 INCKEASE NUMBER UF UNIT=CONVERSION OF TIONS 
SULKCE UATE bHel205 CORRECT TWOeHMANCH CUMPUTED GU TO 
ALSO MAKE AMBIENT=PRESSURE@THRUST=CORRECTION 
DEPEND UN SONIC SPEED aS WELL aS DENSITY 
SUUMCE UATE 68.0828 CORRECT NOZZLEMAMEA SCALING 
SuuRce VATE 6420626 ZEKO OUT NEXT=UNSPECIFIED MANEUVER PRaSE 
SuUKCE VATE 682031} CLEAN UP DRAG PARAMETERS 
SOURCE LaTe 6400301 CHANGE SPELLING OF UNIT =NM] = 
SOURCE UATE BH el2uy9 CONVERT TY CLC 6600 
SOURCE DATE 6700912 FIX INPUT FOR ) = 3 
SOURCE LSTE 6760826 USE CP(7) FOR LAUNCH TIME 
SOURCE OTE 67460726 CALL OLDATA ON KEADING SPECIAL DATA=CARD 
SOuRrt LATE 9720707 ALLUW OPTION M285e ALSU MANEUVER IS IN 
DEGREES IF -FLAG= 1S NON@GLANK 
Suumct Uale 6700263 FOR M)e}» STORE JNPUT MASS IN CiNee) 
SuuRKce Ustt 6601222 REARRANGE OUTPUT» ALLO MANY STAGES MNVR 
SUUKCE UATE 96.COO! ALLOW DWAG PROFILE TU RE JNPUT 


WEAUSs FHINTSe AND SCALES FLIGH! PARAMETERS 


TITLE = MULLERKTTH TITLE TO PRECELE OUTPUT 

“OnE = CUPTROL wURL OF 3 DIGITS Mj oM2eM3 wHICH GOVERN READY 
IN UF THRKUST@MASS PARAMETERS: AERODYNAMIC PARA 
PETERSe AND MANEUVER SCHEDULE®s KESPECTIVELYs AS 


FOLLOWS & 
value PURAMETERS WEAU 
ae 2 ee SO F&O OF BF AES SEARED SSAEHRUDSARABDSR*ENAOMOABETHOGVEDREMOETS 
w[sd © (8GiWE) 
1 MASS, POUNDS (FLIGHT 1S UNPHOPELLEL) 
2 VACUUM TeRUuST¢ POUNDS I[CUNSTANTSTHRUST PROPULSION IS 


UPDINARILY USEGe ALTERNATIVELYs CONSTANT= 
ACCELERKATIUN PHUPULSTON MAY BE HAL BY MANAING 
VetS CARO IN COLS 1°30 INITIAL VACULm THkust 
ACCotGSs AND ENTENTNG APPHOMRIATE VALUES) 

hOCCLE EAIT AWEAs IN SU 

INTTLAL ASS» PYUNDS 

FINAL MASS» POUNDS 
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EGUIVALENCE (AEGe TEU) 
EQUIVALENCE (MAeKM) 

LOGICAL NaMto 


Kakstn > LF EINPCLT.0) RARS) HO 
INPTSSTARS (INP) 

NO 22 INPLTSLeINPTS 

WEAD(Sel) IMAGE 

FORMAT (Tall) 

IF (CHERKFIL (52) KETURI 

IF (IMAGE (1) eNE eHMOLD DATA) GO TO 6 
CALL GLDaTal3nIN}) 

WE TURAN 


® IF(IPetQ.0) GO TO 30 


CECONE (OuseSe IMAGE) NAMyNAM2 


S FORMAT (S0X0A30A7) 


0 2 nel NhV¥ 

TF (NOmMeEQeoMAMVAL (2 0K)) GO TO © 
CUNT TP UE 

DECONE {70039 IMAGE » KoxXMeayD 
NAME IQ 

NAm2e) 

NAMED) # eFALSEo 

Go Tn 10 


& MASNAMVAL (RUNITS¢20n) 


NAME(S. TALES 
AUS0. 
DECUNE (50030 1TMAGE) X 


3 FORMAT (O0ACIF1000) 


agnr 


TF UNM eNE (te) ABACKM 

IF (xLeNEoDe) ABA/XD 

IF (LINE e¢} eGEeMAXLN) CALL MEAD (4Q) 

LIME SLINE 4) 

IF (xeNE eA eORe NAMED) GO TO 19 

WHITE (AGTH) CIMAGE (K) eK] 0%) oh eKAn 
Gu TO 


21 
18 FORMAT(L2OXAGAL OOF LOebeAle 


ADOATCTA Lm(G16e7015mM AFTER SCALING) ) 
WHITE (6018) CIMAGE (A) oA8] 94) eR eKAW eo NAMONAMZ 0A 

Theatre eGTe0) GU TO 22 

te€qea 

asaEo 

OatAcTAPLT) Sa 

RE TUWA 


TEQeTMaGE(5) §$ aBate 
CaLiL COUNOUT()) | 

WHITE (Go3l) IMAGE 
FORMAT (122601004234010) 
6L Tn 22 


Eni IN) 


69.nTS9 
69.0709 
66.0527 
68.0527 
68.0827 
68,0827 
68.0527 
68.0527 
66.0527 
66,0527 
68605c7 
69.0709 
69.8769 
69.0709 
69.0709 
00/15/49 


68.057 
69.0709 


68.3205 
68.0209 
976 


66.0209 
68.9209 
67.0726 
67.0726 
67.0726 
66,0507 
6920709 
69.07U9 
68.0527 
68.0527 
67.0712 


69.0709 
69.n709 
o7.071k 
67.0711 


69.0709 
67.0711 


67eN711 
67.0711 
67.071) 
930 
9@) 
962 
943 
69.0709 
66.1205 
9e7 
66.1205 


69.0709 
69.0709 
68.209 


991 
992 


68.0507 
69.0709 
68.0S07 
69.0220 
68.0807 


SUBROUTING INA(DOATAONTN) 

SOURCE OATE 68.1208 PRINT ASTERISK If DATAsvALUE CONVERTED TO INTEGER 
QOURCE DATE 60,0209 CONVERT TO CDC 0409 

@OURCE DATE 67-0726 CALL OLDATA ON MEADING SPECIAL DATacCand 

SOURCE OATE 67.0713 ORaNN new COOE 


READS NIN VALUES FROM IN1L°TYPE CaRMOS, ANDO STO@ES THEM IN THE Data 
aRRAy, THE OLFFERENCE BEING THAT IND READS 3 VALUES FROM 
COL 41070 OF GACH Cand, 


aananaaanen 


000005 MIMENSTON DATAININ) 0 10(6) oVAL¢3) 
000005 POUT VALENCE (AE Gs EQ) 


000008 waReyn & ITFCMINLY,0) KARE] NO 06,1208 
000031 MAX @ TABS(NIN) 
000013 MK @ (MAKOD}/9 
o00ol7 uad 
000020 NO © MXKa) MK 
o0eez) MEAN (Sel) Deval 
000030 PORMAT (441093F10.,0) 68,0200 
00eese TPH (CHEKFIL (59) RETURN 08,0200 
000034 TPIZOCI) NE SBHOLO DATA) 60 To 6 67,0726 
000040 CALL OLOATA(3HINS) 67,0726 
000043 RETURN 67,0726 
060042 CALL COUN OUT(1) 68,0209 
WRITE (608) LO, VAL C1) KAR 68.1206 
PORMAT (1292016 0F14.6,A1) 68,2208 
jr ime]. Of,MAX) 60 TO § 
CALL GOUN OUT«(1) 68,0209 
WRITE (G03) VAL (2) oKAR 68,1208 
PORMAT(F66.4eA3) 68.1205 
17 (M02,08,maK) @0 TO § 
CALL COUN OUT«(}) 68,0200 
WAIVE (03) VAL (3) okAR 60,1206 
nO © Rel. 
4 @ Mel 
TP(M, OT MAX) RETURN 
IFI(NIN,OT.0) 60 TO 6 
TEQwvaL(ny 8 yal (m)eage 
000140 HATAc(m) © VAL(K) 


000350 RETURN 
00015) eNO NI 


omen 


a gaa semilinear 


: INCOL 
, © 
: 
& 2 
‘ 
é SURRCUTINE INCOLETITILE e SETUP NAMES eUATaeL INES? 113) 
: C Soumre VATE 6900709 WIVEN Fo3 ANU Feo FURMATHFIELDS 
; . C SLuNrE UAlL 68005U7 PROVIDE FOR UPENSENG KEALING 
2 . . Scunct DaTk eH.0209 CONVERT TL CLC e800 
é c Sucmet LATE @7elU27 WEPATK CHANGES OF 6720726 
% Cc Scurce OATE |7eC?Z6 CALL ULDATA ON WEADING SPECIAL DaATaeCako 
C SOumct ualt 67eMa0o COMWECT- D FMTS ARKAYs KOLUMS COUNTER 
c Suukce DATE 65e0i06 
c 1132 
. c ~EAUS AND PRIWTS COLUMNS OF DATA FROM CARUS 1333 
1136 
) é AMGUMENTS ARE DEFINEC In SURMOUTINES OUTCUL AND OUTSET 1235 
C 1136 
c AULICGL PRINTS LINES OF 1 TO 10 NUMBERS IN COLUMNS, af Tn 1137 
c AEALTINGS aN TITLeESe INCOL USES OUTCUL FUR THIS PURPOSE 113A 
c oe TER FIRST HEALING THE LINES FROM CakDSe IN DOING TeISe 1139 
Cc 1T CCASIUERS CakUuS TO Se UIVILEU INTO 7 1U*CHAKACTEw COLUMNS, 1140 
Cc SO IF & LINE CONTAINS 6 Om MOKE COLUMNS TwO CARDS ARE KEQUIREU 116) 
Cc Pew LCINEs INCOL CETFRMINES FORMATS FOR EACH CARL COLUMN FROM lle2 
c FAAMIAATION OF Tre FORMATS IN SETUP TU BE USED BY OLTCOL IN 113 
¢ PKINTING Tet INPUT Dala&e lles 
c 1165 
c TNPUT DATA IN & LINE IS LOCATED ON CakUSS) aS FULLOWS - 1146 
1Lle? 
TAR } . COLUMNS 1-10 ° WORD }) 1] 468 
¢c 11°20 = WORL 2 1169 
¢ aND SO ON FUR BURDS 3 = 7 1180 
Cc caRD 2 © CCLUMAS 1-40 = {NOT wE~AD) llsi 
¢ 41°50 = WwOKO @ 11se 
f @ ¢ 51-60 = BOND Y 1183 
‘, c 61°70 | BOND 10 1156 
r c 1155 
i c Tr At NeCHARACTEW NAME IS REQUIHED FOR TmE LINE® {7S N CHaRe 
. c ACTERS APPEAR IN COLSe 1eN OF CARD 1] AletD THE wORUS ARE MOVED N 
es c rakO COOLS. TU THE RIGMT IN CONSEVUENCEe 
: c 1159 
Cc FORMATS FOR DATA CULUMNS IN EACH CARD ARE DETERMINED FRUM 1160 
c FURMATS IN SETUP OK TRE CODES SELECTING TrEM IN OUTSET 1iel 
C AS FOLLOWS © 1ie2 
Cc hle3 2 
Cc rae 4 GUTPUT FORMAT INPUT FORMAT lee 
Cc oor woecnewncscr otwesccoenes 1165 
c i El2e4 Elveu 
¢ ? Elle2 Elved 
¢ a Flleo FlO6c0 OYeNTLD 
Cc 4 Flie3 Floey 69-0769 
c & Fi06d Flueod 
( 6 18 Fluev ~y 
c ’ 119 Fluceu 
. c Ry O12 Olo 1173 
¢c 9 alg Alo oH.n2(9 
¢ , AGea2 Abea2 
c 1176 
c CECIMAL PUINTS PUNCHED IN THE CARLS OVERRIDE, OF COUmSEs 1177 
c TROSE IN TRE FORMATS. IF THE OUTPUT FURMAT IS INTEGEHs THE 1178 
g c IAPUT FLOATING PUINT VALUE IS FIXED ACCORLINGLYe 1179 
¢ 1Leo 
ce . ¢c th Tre FIRST wORD OF TITLE IS AN ASTERISK FOLLUMEL BY BLANKSe isi 
. e ¢ GUTPUT TS OMTITEN. ws2 
c 1103 
c whITTEN 1/8765 1164 
: € 1145 
€ CCOMPKG,FHHATS 69-0709 
4 | 000019 COMMON /ERMATS/S PEORMS (2913) 
j 0 ATA (PFUNMS 8 
; i o00ul (aReyeies el2nG,l elk ol QMFL L601 ’ 69.0709 
P| 2 1PKF1] 03014 #12HF 10 6u92X oL2HIOs4er ‘ 69-0709 
j 3 1amtioeax el2nol2s vo L2HAALO OX ’ 68.0209 
4 @ & LeMiAABeA2elA POL2HAXZABeARS2K SLAHIAZABSAGSIA 2 
S LomaxSabeak H piaiee 
OMPKG OF MMATS . 
| c ae 09715769 


ake Ss) 2) ee 


mnie. 


Pos nneag 7 eS 


COOuIS 


0000.0 


o000)° 
o00vld 


oo00nl0 
000011 
o00n1) 
000014 
o00u1s 
o00n20 
000022 
ononern 
000027 
000u36 
000035 
000036 
000000 
000062 
000466 
000067 
O00use 
000052 
oo00rss 
000087 
000086 
000n7sS 
HAVO??? 
000101 
0n01¢62 
000103 
000107 
O00T1} 
00113 
0091315 
000116 
000123 
000123 
000127 
000) 36 
0001 %6 
000137 
000137 
00914) 
0001467 
000153 
000154 
0001480 
0001646 
00026 
000216 


000266 
000266 
000257 
000284 
000256 
000260 
0002762 
000265 
000267 
O0027TH 
000277 
000303 
000303 
000306 
000305 


4-74 


1? 


le 


UTMENSICN MTS (16) SCAMDFM (1S) pSE TUM (G2) eNAMES (502) sUATA(2) 6 


e TNTOR (YG) oP aGE (8) oTmaG2(e? O9eNLGY 
LATA COPMTS UT otmdbelede 21OmMEL De Oo) eS1OMFLOe0e ebHOl0s 8 

1 OFaiie APA AZ ee BM 2AB Abs eBMJAbLAO, »OMSAE, Osean 6 68.0209 

2 felele)) 


CATR (BLANBF®I Mm 1s (aSTRSABI Me) 
EUUTVALEACE (abwo IEG) 


Laval Ines 6620507 
IF ULIMESsE G00) GC TO 14 65.0209 
CaLhL AMET ColN eat NTGR) 67.0406 
Corurm (Li SER TS (LS) 

Nae NeRULLMSsO 6721027 
KEDCOL ENR AWD 3] O7e1Ge! 
vem eCOLBlely 1200 


OU ¢ WPRINTZ1 613 

TF (SETUP (2@MCCL) -EU.PFURMS() oMPHIMI)) FO TU « 

CONTIALE 1203 
6c TO A 12fe 
KHUCOLBKREUCCL ©] 

IF (MPRIATLT 11) GO TO 4% 

NAMEN SPPRINT on 

CarUr rh (RmUCCL) BETS (MPHING) 

bU S TR) LeMPHWINT 

RRECOLMARLCCL ©} 

CARDFEM(RRLCOL) SOLANA 


GO TA A 1216 
KULUeSBACLUMS*) © MPRINT/LO O7enade 
Th (PRINT Ewe Che MPM INT cE Ge?) INTGOM(KOLUMS) &) 1206 


{F (anOCCLehtev) GO TO 7 
CORUFR (KwUCCL BPM IS (Le) 


RPO COLeRRE CLE OL 1219 
NAAR SQ 67.1027 
CARUEM IR 6LCCL) SFR 1S (MPRINT) 1220 
Cor Trae 122) 
CARDE (RnUCCL ©] ) SF MTS (36) 

LYPESBTAnS (IL Pk S) 1223 
LU 18 LYNE SLeLYNES 1226 
HEAU(Sel) ThAGE 67.07¢6 
FORMAT (HALO) 66.n26Y 
TFC CHERFIL (50) mE TURN 66 0N2C9 
Tf (IMaGe (1) eEUsBRMEND OATA) GO TO 20 

TF (TAGE (1) eNEsBROLD DATA) GO TO }} 67.0726 
hore Siteol) e1MaGE (1) OT eNT2H 
we TOW, 67.07 CR 
TP (NKARDCEUel) GO TU 33 

WEAUIS el) TAG? 6761027 
TF (CmeRP IL (59) KE TUN 68.02UY 
T2mLvnE@KCLUMS 

TyslaenCLumse) 1226 
IF {th OPELeNE eG) GU TO 12 1227 
DECOE (on eCARDFRolMAGE? (CaTatT)elellel2) 67 eNTe6 
GG Tt 148 1269 
DECUTE Cai CARO Ms IMAGE) (NAMES (JeLYNE) eJeleNartu) » 

° (VATACT oleliel2) 6TeN726 
KOLUMEn 1231 
bu 16 Tallel2 1232 
RCL Ue KUL UMO] 1233 


IF (1> TGMCRULUM) e600) GU TO 16 
TEGSNATACT) 
CATAt(I) BAEG 


CORT IAUE 1235 
CONTTAUE 1236 
Th (TLE NES ASTHSK) CALL OUTCOLI TS ILe oSETUP Name SeLaATacL Sav) 68eNSC? 
RETUPR 1236 
LSAVeISIGNILYNE@} OL INES) 66,0507? 
LINESSLYNE =) 68.0507 
TF (LSAVenNE oO? GO TO 19 68.0807 
RE TURK 68.057 
ENG INCOL 


— 


p 
y 


SURMOLTINE INUEC (ahGeMm) 

Scumct VATE @ye0TU9 REMOVE ENCODE STATEMENTS 

SULmCE LATE 690013) aALLU® PRINTOUT TO BE SUFPHESSEL BY CALL INLOU10) 
Suudrk t-Te #60209 CONVERT Yu COC 6000 

Suumct valf obevocl 


FUCRMAT@FREE [MPLT RUUTINE 


FrTmeres InDEC RETURNS FLOATING*PUINT NUMDERS 
Intnl KE TURNS FIKEDOPOINT NUMBERS 

MEALS DATA CamUS UNTIL MM NUMBERS ARE FUUNC® 
WHICH ARE STOWEL IN AkGe 
SCANS OaTa CAWOS UNTIL A LEGAL NUMBEROFIELO 
HEGINSe Tmt CONTINUES SCANNING UNTIL 
THAT NUMBEK@6 1ELD TEMPINATES ON AN ILLEGAL 
CHARACTER UM A BLANK © AND REPEATS mm Times 


IAVRAINANN ANANDA BEAN 


an INTEGER FLELU FULL OwEY, IMMEQITATELY AY 
A LEFT PAwENTMESIS 15 INTERPHE TED aS @ 
BEPETITION FacTOR (Set SamMPLE DaTa=CaRDd) 


SAMPLE DATA CAKO woe 


KeVoZ COOKDE2e30FS GIMEL TO HaDAR) 060000 21100) 
NUMBERS eeecee @ eeeee eve 


WRONG DATA CahD coe 


WeVeZ COORLB2c3@ ES LikEr TU WAL 2) 600011000100 
coe PROBABLE ERROKS © ° e * 


RETURNS MM WONLS UF BCL Lala 
STARTING al FIST NON@bLANK ®ORL 
@ ORO 3S LEEMEL TU END ON & BLANK: & COMMAy 
OM THE YENTH CHARACTER 


“ese TS SO Treat neat CaLe UF INDEC OR ININT 
WILL START ON New Cart 


OTR eos NORMALLY LATA CARDS ahE KEAD FROM 
5 AND PRINTEL Of 6 
ENTRY IPRLIK CeanGES INPUT UNST TO mm 
EN Try TAL CL CranGES OUTPUT UNIT TO mm 


¢ 
& 
¢ 
c 
c 
C 
c 
¢ 
c 
¢ 
¢ 
- 
é 
C 
- 
c 
¢ 
¢ 
c 
¢ 
Cc 
c 
c 
is 
¢ 
¢ 


000005 GIMENSTON RAmRU LE) e ARG (MM) 69-0709 
0000u8 CaTarmOUNTSO ela] BO) ef 7Pe SioellOs 6) e (MULTSO) & (NCPWSI 0) e (KaAR@O) 69eNT0G 
ooonus EurcTval€NCe (AEG IEW) 09.0709 


0000n% Tr Toe kee s GC TO 14 

O00%UK EStTry [hTT 

o00n0ts INTUOFH®) 

O0Cl 16 1G Ctl ae MBL oe MEY 

000020 Te (ADUNT eC GTeR1) GU TO 20 

000023 e eau TPe)) Rand 

000031 FuRmal(84)¢) 69.0709 
000031 Te (CME RP TL CTR) RE TURE 69.0131 
000035 IF C1%eE em) CALL COUNQUT EO)? 69.013) 
000063 TROD Ge) wr ITELTOe3) RARD 69.013: 
000054 FOwM*T (Sanal fy} 69.0709 
o00css AUUNT SHO 

000058 K126 


0000546 Ni Poel COTAMY INENFLSSNFE SO 

000063 Klente) > It (K10GT28N) GU TO 12 

00007? TE (I n1GtreL leu) GU 10 $0 

000073 CALL MCMARI(K] pRARU gNCPWweKAMe 1) 69.0709 
o00n7T IF (Ramek Gel mo) ARI NOY 

000102 IF (RameECodRe} NDOTS) 

o001uSs TF (KahbeGtelHO cANDe KRaReLeoIkG) NOIG9)] 68200209 
000117 It (NW INCH DOTONDIGeLE SO) GO TO 22 

000126 R2eK} 

000125 MINE WSTONMIA 66.0209 


2 OREO eT ANTUNES 


OONIT] 
060375 
0004u6 
000407 
00061) 
000415 
O00e2u 
0094633 
00n435 
000037 
O00N4eu 
O00Gen 
000650 


00045? 
000s) 
000463 
000472 
O00H74 
000503 
0005%6 
000505 


000504 


000126 
000333 
000136 
000147 
0001S 
000151 
090153 
000}6) 
000163 
000171 
000173 
000203 
000205 


0002146 
000216 
006225 
000231 
000238 
000236 
000237 
000247 
000260 
000264 
000276 
000304 
000316 
000327 
00033) 
000336 
000337 
000340 
000341 
000343 
000345 
000346 
000351 
000383 


000383 
000362 
000363 
0003466 
000346 


4-76 


AA 


S¢ 


se 


Se 


36 


36 


4? 


Ala 


“0 


46 
“a 


4“? 


09 


CALL PCNA CK RAR NCPmORARG]} 

TP (Rab eeu e im Ome RAW EGLIRe} GO TO 09 
00 SP Ke2eNCPa 

Kean} on=] 

JF (RP 6GTeRC) GU Tu Se 

CALL MCrAm(K2 RAHUL eNCPweKARe}) 

IF (RAW SECeLH Oe KAR EQel Ke} GO TU S@ 
KR2eKk} eNCPa 

K2anpe} 

TEQ8)h 

CALL MCHAR(IK] hARGel eo TEQoR2en] 1) 

K 18K? 

6C 1. 06 


ERTHY TNL IN 

TPemm $ GL TO 90 
ENTRY JNLOU 

[Gmamm § GO TO GU 
ENTRY INAEte 
KUUNTaN 

MULTaf 

Age TURN 


ENn  JADEC 


KeeKkre) s 1F (K2eGTedn) GO TO 3e 

CALL MCrram IAS eRAwD eNCPReKAne)) 

LF (RAR LT eI RD COMeKAN.GTOINGO) GO TU 26 
KLIGahDIGe) 

GO To 26 

IF (ROW eNE CIT Re? GO TU 28 

NUOTSNOOTe] s IF INDOT¢GT.1) 3e02e 

IF (RAR NEC] RE} GO TO 30 

NEESPEEC] $ TF (NEE SGTol) 30026 

Tf (KAR NES} HS) GO TO 32 

NPL SeAPLS*) s IF (INPLS®°NMINOMINE RSToGTe}) 3e020 
IF (MAW eNESIR}) GO TO 36 

NMI NSNMING) 8 IF (NMINOMINFE RST.GT.1) 34926 


IF (NNIGeLEeG) GO TO 20 

TF (MULT eGTe0 eORe KAKeNEIRE } GO TU 36 

TF (NDUTONEE °NPLSONMINGEQ20) MUL T@~) 

K2an2~} $ KaaK2en}eo) 

FeTaar( t ) 

IF URteEUel) GLU TO 37 

Le(Kre1)710 & CALL MCHARINCPwelolkusceFmte]) 

Laklel-lNeL % CALL MCHARINCPweLolkGedeFmTe)) 
CALL MCHAMINCPwWs. IRAs GoFMT9)? 

LeKR/10 % CALL MCHARINCPweLolROseeFMT5)) 

Lexkeln®L S$ CALL MCHARINCPWeLeIROeTeFMTe}) 

DECONE (K2eF MTs KARO) AQ 

TF CINTGER SGT 9) eOK, MULT LT 20) TEQeAtY 

a) eK? 

IF (MULT) 38206060 

KMUL@K18K1¢1 

MULT@IEQe) 

Ge TH 20 

IF (RaW eNEC}M) ) GO TO 46 

KI} eKmUL 

MUL TamULTe} 

AWG (~) BAEU 

CONTINUE 

GO Te 9 


EuTRy JNALF 

INTGERS=) 

GG To 10 

Ky} akyeo} 

TF (K}eGT.HC) GO TU 12 


692709 
69.0709 
O9-n7TU9 


69.NTUY 
69.0709 
69.0709 


65.0709 
69.0709 


OF-607( 4 
686.0269 


Ow oNZLY 


OB.0209 


69.0709 
65.0709 
649-07CY 
69.0709 
69.0709 
65.0709 
O900709 
69.0769 


69.0709 
69.0709 


FUNCTION TNOEQU (LENToe MARAV, MATCH) 

SOURCE DATE 6620209 CONVERT TO COC 6400 - @RAND NEW CODE 
FOR ENTRY INDEQUs RETURNS INDEA I SUCH THAT maRAY( 1) MATCH 

FOR ENTRY INDMAX, RETURNS INDEX I SUCM TmaT maRayil) IS LaRGEST 
FOR ENTRY INOMIN, RETURNS INDEX 1 SUCH THAT mARAV(I) IS SMALLEST 


ee00es DIMENSION MAKAYILENTH) 


000005 00 12 Lal.sLENTH 
IF (MARAYIL) cEQ.mATCH) GO TO 16 
CONT INUE 
20 

eeeols INDE QUSL 

000016 RETURN 


eoeols ENTRY INDMIN 
LLsO0 8 IFILENTHO]1) 26021923 
LL#] 8 60 10 26 
My AL eMaRAy()) 
Lai 
DO 26 LEZcLENTH 
IF (MARAYI(L) eGEoMVAL) GO TO 24 
MVALSMARAY(L) 
LLeL 
CONTINUE 
TNOEQUSLLE 
RETURN 


ENTRY [NOMAA 
LL#O 8 IFSLENTH©1) 36033033 
000071 LL9h $ GO TO 36 
000073 MyALeMARAY(}) 
oeeo7e tL*) 
000077 00 36 Le2eLENTH 
00010) TF(MARAYV(L) eLEoM¥AL) GO TO 34 
000106 MVAL@MARAYV IL) 
0¢e0110 Lest 
000111 CONTINUE | 
000114 INDE OUSLL 
000116 RETURN 
000116 ENO INDEQU 


elma. 


SUBROUTINE INDOCTITLE sLOOPS emaSE VALUES) 

SOURCE DATE 68.0717 FIR LOOPING IF DOCON (a. » © BLANK 
SOURCE LATE 68.6209 CONVERT TO COC e000 

SOURCE DATE O700726 CALL OLOATA ON READING SPECTAL DATACARD 
SOURCE DATE 6700726 BRAND nES CODE 

SOURCE DATE 6600601 CALL OUTSET wiTm 300 Im FIAST ARG 


INPUTS, OUTPUTS, AND INCREMENTS UP To 1¢ GENERALIZED 00 LOOPS 


TITLE WOLLEATIM TITLE POR CUTOUT 

LooPs NUMBER OF LOOPS TO RE INCREMENTED 

waS€ NUMBER OF TRIAL IN PROGRESS 

VALUES LIST OF VALUES oF INCREMENTED VaRtaaces 


FOR KaSEwele INDO READS ANO PRINTS FOR Each LOOP a CARD 
SPECIFYING LOOP INCREMENTATION © 


COLUMNS coment 

@eececeee @eeceeare 

1730 WOLLERI TH LOENTIFICATION OF (OOF 
VARTAGLE 

31-40 F10.0 MINIMUM VALUE ASSUMED OY VARIABLE 

#1°S0 1000 MARIMUM VALUE ASSUMED BY VARIABLE 

$1°60 71000 VALUE SELECTOR FOR vaRIaaLE 

61°70 F1000 SELECTION METHOD CODE 


SELECTION METHODS ARE DEFINED AS FOLLOWS @ 


“E TROD VALUE SELECTION 
cooEe SELECTOR mETHOO 


0 OR 1 INCREMENT LINEAR INCREMENTATION 


ERAMPLES cee 

FOR LOOP FROM 1 TO 3 IN STEPS OF 15 
PUNCH le Je «3 le 

FOR LOOP FROM 3 TO 1 IN STEPS OF 2S 
PUNCH le 36 e658 Le 


STEeS/ LOGARITHMIC INCREMENTATION 
OECADE 


ERAMPLES woo 
FOR LOOP FRom 3 TO 300 at & STEPS/ODECADE 
PUNCH Be 3006 @e 2e 
FOR LOOP From 300 TO 3 AT & STEPS/DECADE 
PUNCH Je 3004 he 2e 


(INCREMENT) RANDOM SELECTION 


IF VALUE SELECTOR IS ZERO, VALUES ae 
TAKEN aT RANNOM BETWEEN MIN AND aX VaLUFS 
EXAMPLE coe 
FOR RANDOM VARI AQLE aa S aNn 
PUNCH Se 6e Oe Jo 
TF VALUE SELECTOR IS NON@ZEROe VALUES ane 
PICKED RANDOMLY FROM THE SET OF OTSCRETE 
VALUES WHICH WOULD OCCUR UNDER METHOD 1 
EXAMPLE coe 
FOR RANDOM INTEGERS RETWEEN 1 AND 10 
PUNCH te 1@e te Je 


& VARTABLE SELECTED BY METHOD 3 1S NOT, STRICTLY SPEAKING, 
IN & LOOP © THE SEQUENCE DOES NOT TEAuINATE ITSELF. auY MUST 
BE STOPPED BY THE phd ad bed OF OTHER (TRUE) LOOPS: OR BY 
OISCONTINUING THE CALLS OF INDO. 


¢ 
. 
¢ 
c 
c 
c 
¢ 
c 
c 
c 
¢ 
c 
c 
¢ 
c 
c 
€ 
¢ 
¢€ 
€ 
c 
€ 
c 
c 
Cc 
¢ 
¢ 
c 
c 
Cc 
c 
c 
€ 
¢ 
€ 
c 
€ 
€ 
c 
¢€ 
c 
€ 
¢ 
c 
¢ 
c 
c 
lo 
if 
c 
¢ 
c 
c 
€ 
€ 
¢ 
c 
c 
€ 
¢ 
¢ 
€ 
c 
¢ 
c 
c 
¢ 


VARIABLES SELECTED SY METHOOS @. 16 OR 2 aRE ORDINARY LOOP 
VARIABLES © THEY PROGRESS FROM One EXTREME VALUE TO THE OTHER. 
16 Tee METHOD SELECTION CODE I8 POSITIVE, THE LOOP VILL BE 
WESTED = THE OUTER“LOOP (EARLIERREAD) VARIABLES WILL NOT 
CHANGE UNTIL THIS CINNER) LOOP Ig SATISFIED: FOR NEGATIVE 
METHOD CODESs TurS RESTRAINT 18 RENOVEDe ALLOWING SinuLTangous 
STEPPING. EXAMPLE coc 

OATAcCAROS PUNCHED (A) Yo 2e 0S te 

(OB) 10. 1000 2e @2e 

RESULTING SEQUENCE (eB) @ (leol@ede(leSeBleb) o(2eel00.) 


PHEN LOOPS ARE NESTED. THE LASTeREAD LOOP IS INNERMOST 


FOR KaSE GREATER THAN OR EQUAL To ZERO, INDO INCREMENTS KaSe 
BY ONEs SELECTS NEW VALUES AS REQUIRED. IF ALL LOOPS aRe THUS 
SET aT FINAL VALUES: INDO RESETS KASE TO 6. 


TYPICAL USE IS aS FOLLOWS ~ 


1 KASEse} 
CALL INDOITITLE cL. 00PS en aS€ eo VALUES) 
{10 READ OATAqCARDS) 
2 CALL INDO(TITLE cL 00PS oKaSE oVaLUES) 
(TO SET UP VALUES) 
IF (maSEcEQ.0) 60 TO 1 
1F (WASE eGE MAX) STOP 
(COMPUTE ® PRINTs ETC) 
GO To 2 


If TITLE IS AN ASTERISK FOLLOWED BY BLANKS» ALL OUTPUT IS 
OELETEDe 


COMPKG.»BASCON 07.0726 
000006 COMMON /BASCON/ 
1 PROGRM, KPaGEe LINEs TOF HAYs RUNe RUNID (66 
2 MSGe FLaGe OATEs MAXTMe MAKPG. MAKLNe 
+ 3 KLASSe KGROUP, KUNITS. KSTINTs KOORD > IFEOr 
@ oNDFILE 


1287 
000006 DIMENSSON VALUES (2) 06 (02) oNAMES (S510) sDOCON (4010) 1268 
000006 OIMENSION VF IRST (10) eMODE (10) 0AA() 0) AB(10) oNTMC LOD eNTIM( 10) 
000000 OaTaA (aSTASKel He) 
000006 OIMENSTON 00C (40) 67.0726 


Ge0e0e IF (MASE) 10002000300 
CALL OUTSET(34G012M MINIMUM =6sleH | =6VALUE 
CALL OUTSET( Se012M MAXIMUM olen VALUE 
CALL OUTSET 64012" VALVE elgm SELECTOR 
CALL OUTSET( POo12M SELECTION o12m METHOD 
CALL INCOLI TITLE oe nAMES eD0C 0 0OPg) 
TF(NOFILE) RETURN 
TF CNAMES (1019 oNEoBHOLD DATA) 6O TO 1¢2 
CALL OLDATA(OnINOO) 
GO TO 106 
CALL KMIT(4@000C e00CON) 67.0726 
CaLl OUTSET (386012" olow 
00 106 Lei eLoors o7.e717 
TF (OOCON (Ooh) 0EGs90-) DOCON ack) S06. O7.0717 
CONT INUE 67.0717 
KaSEs0 1207 
eoeo7e RETURN 1206 


eeeoT7s ASSIGN 236 TO KUMBAK 

DO 236 LPs1sl00PS 
OVSQOCON( 304”) 
avedes (Ov) 
ODeDOCON (2 oh P) —D0CON (1 LP) 
vveQOCOm (1 oF) 
TF (OVeLTe@ed V¥EQOCONIZeLP) 
WPIRST (LP avy 
MODE (LP) 2aBS (DOCON (4 LP)? 
1F (MODE (LP) cEGcO) MODE (LP) 82 
WeMOOE (LP) 
@O TO (21ee2Be0Z30) om 

210 aatLPreov 
SO(LP)el. 
NTMILP) OF DIV (OO/AV) ©6909 
@0 To 222 

220 AAILPI OO. 


000250 
e026) 


000263 
c002e4 
000265 
000267 
00027) 
000274 
900275 
000310 
ee0312 
600316 
000320 
000322 
o0e326 
0327 
e327 
0335 
000340 


000340 
000350 


900353 

000356 

900366 
e 


000033 
000003 


C00ee4 


222 
230 


232 
234 
236 


300 


302 
303 


304 
306 


000 


se2 


900 
902 
906 


OB(LP)010.%°FOIV( 1./0¥) 

NTM (LP) OAL 062 0 (FOE Y (OOCON(ZeLP) /D0CON (1 0LP))) Cae .089 
VALUES (LP avy 

GO TO 236 

IF (OVeNEsOe) 60 TO 232 

AAILP) SDOCON (1 LP) 

86(LP) sD0 

60 To 236 

AatLP rev 

OR (LP) SFOTV (00/av) 

NTM(LP) 80 

@0 To 500 

NTIM (LP) OTSIGN (NTH (LP) sOOCON(a0L PF)? 
GO To 900 


ASSIGN 300 TO KUMBAK] 

ASSIGN 304 TO KUM@AK 

00 306 LPa1 LOOPS 

LLeLPos 

IF(LLGTeLOOPS) GO TO 303 

00 302 LeLteLOoPs 

TF INTHIL) cGTeO cANDe DOCONLOst )c6Tede) GO TO 306 
CONTINUE 

OVEQOCON( 30LP) 

IF (MODE (LP) 06063) GO TO S00 

IF (NTMILP) oME00) GO TO 000 
VALUES (LP) OVFIRST (LP) 

NTM(LP) ONTIMILE) 

60 TO 306 

NTM(LP) OISIGNC LABS (NTMILP) Do) oN TMCLP)) 
CONT TINUE 

60 To 900 


VALUES (LP) OVALVES (LP) COB (LP) CAAILPD 
60 To KUM@~AK] 


IF(OVeNECOe) GO TO S02 

VALUES (L®) CRANE (90) CBB (LP) CAR (LP) 

60 To KUMBAK 

VALUES (LP) COVOAINT (RANF (OQ) OAINTIBBILP)OLoI}¢ AACLP) 
@0 To KUMBAK 


DO 902 LPa) sL00PS 

IF (NTMILP) o6T20) GO TO 900 

CONT INUE 

KaS€ac} 

KASEOKASES) 

a ate cel CALL OUTCOL (TITLE oe NAMES» VALUES cL COPS) 
€ TURN 


68.0200 
o8.ezee 


6a .0209 


4-81 


0000uU6 


000006 


000004 
000006 
000004 


Sa pl ean IRA ae 5 % 


Q0000A 


000604 


SURMOCUTINE [Nav 
Source LATE 6920709 
SOURCE UATE 0401205 
SOLACE DATE 6829720 


Inmv JTAPLTS a MOLE AND 
CaRD FOMMAT «ee 
Cat $1260 VALLE 


UNIT TeNAME PUNCHED 
O&G 
FT 
Kw 
NMI 
KFT 
M 


aAaANAN DANN ANAANAAANANAANANAANA 


COMPKG.RASCON 
COMMON /HASCON/ 


1 PROGHMs KPAGE® 
2 MSGe FLAG» 
3 KLASSe KGHOUP» 
& oNnb ILE 
CCOMPKG+BASCON 
{44 
CCOMPRKGeCONCON 
COMMON /CONCON/ 
1 Pie SRUe 
2 GacCe GCONe 
CCOMPRKGeCONCON 


¢ 
CCOMPKG se CONVRT 


(TITLE eMOUE oe VALUE) 


WEMOVE CALL UF AMIXER 


aQD COMMON /6ASCON/S 


BRAND NEW COLE 


& VALUE FROM a UaTa CakO 


TITLE IS & BCD CHARACTER STRING ENDING UN A 04SHA 
COL lee) ANY COMMENTS 
CoOL #1250 MOLEs PUNCHED IN FeFORMaT 


COL 61263 UNITS OF VALUEs PUNCHEL IN A@FORMAT 
VALUF WILL BE SCALED TO INTERNAL UNITS AS SHOWN BELOR coe 


INTERNAL (MET.) (ENGL) 
RADIANS RADIANS 
METERS FEET 
METERS FEET 
METERS FEET 
METERS FEET 
METERS FEET 


LINE®s TOFDAYs RUN? 
Oates MAATMs MaXPGe 
KUNT TS» KSTINT, KOURD, 


SLVe SMP Skre 
wROOYe RRUZROe TwOPT >» 


COMMAN /CULVRT/S NNVS NAMVAL (391)) 


NATA (NNv@))) 
DATA (NAPVAL & 


1 BnVEGe 2(601745329252) 9 
2 2ne To 0.304HU065 1.0% 

3 JHPSF, 4.88024, 160% 
ry 2HKMs 10000s 3280¢R33e 
5 JHNMIs = 19539025e 608Ne3e 
6 2nL8e 0.453592 109 

? 1G, 9.60665— 32¢17398¢ 
) 3xKFTs 304¢8006* 1000e0 
9 Lae leds 30280833 + 
0 BHMRA, 212001) 9 

1 BHSECe 2lled 3} 

CCOMPRG se CONVAT 
DIMENSTON TEXT (a) 


DaATa (ANDe]Pr ) 


RUNID (6) 6 
MAALNe 
IF EOF 


RBUDYe 
HaFPY 


09715769 


09715769 
69.0704 
6920709 
6920799 
68 efS27 
68.027 
68.0527 
66.0527 
68.0527 
6820527 
68,0827 
68.0527 
68.0527 
69.0709 
69.0709 
69.0709 
6900709 

09715769 
69.0709 


2p needa. 


Sse 6 ee Oe 


¥ 
a 
ag 


fire ht iT 


pics ae 


See ek HE wet eet ee = 


a 
) 
3 


000004 


d000U5 
000010 
000024 
d0un26 
000u32 
0000346 
00004) 
000083 
0004s 
000065 
00005) 
0000S? 
000054 
000055 
000001 
000162 
ooolV0e 


000)03 
000123 
000)26 


OU 


lb 


EQUIVALENCE (ARGS TEU) 


CALL TITLEM (ele TITLE) 

REAL (Sel) TEATSODE sVALUSNAMeNaM2 
FORMAT (GAL U02F 10eA30AT) 

Th (CHEKF IL (S)) BETURN 

DU 1a NELeNNV 

IF (NAMeEGeNAMVAL (10ND) GO TO 12 

CONTIMUE 

VALUFeVALU 

GU Tn 16 

Te 8e-aMVal (RUNITS¢20N) 

VALUF SVALU®ALY 

MODE wODE 

CALL COUNUUT (1) 

IF (VALUReEGeVALL) GU TU 18 

WRITE (602) TEATOMUDE oe VALUSNAMeNAM2 9 ANU 0 VALUE 
FORMAT (LSAGAL001A IKE 100303XAZ AT CASS (OG130608 AFTER SCALING) ®) 
RETURA 


WHITE 1622) TEATePUDE s VALUSNAMpNAM2 
RETURN 
ENO INMV 


69.0709 


69.0709 
69.0709 


696709 


69-n709 
09.0709 


69.0709 
69.0709 


69.0709 


te ect < ¥ a Boe sninsiniibtiese - 


a 
JACOBI 
GUBROUTINE JACOBI (AeMeNeACCUs Ven? 
: e @OURCE OATE 60.1205 MAKE 17 WOR FOR DIAGONAL MATRICESs TOO z 
e SOURCE DATE 68.0200 CONVERT TO COC 6600 3: 
| c QOURCE OATE 6e.07EL BRAND NEW CODE 
cuacon; SUAROUTINE FOR Qndodn Poi 06s mosesees 
¢c i 
| ¢ a TS A REAL SYMMETRIC MATRIX 
| ¢ a JS & OUMMY TF T¥O0, IF TVale 18 LOCATION WHERE ELSENVECTOR 
- 7 ¢ MATRIX 15 CREATEDs AWD IF BVez Tg A MATAIN WHTCM WILL BE 
- > C PRE“UULTIPLIED By THE EIGENVECTOR MATRIX, 
| ¢c wy IS StZe oF a ANd 6 
c ACCU TS A PRECISION VALUE. IF LY 1efolGs WILL FORCE TO 1.810, 
=< ¢ TV IS A SwITCu WHSCH MAY TAKE ON VALUES 0, 1. OR 2. SEE @ 
a ¢ PARAMETER DESCRIPTION ABOVE. 
; c wR Ig AN ITERATION COUNTER, 
a e C { 
q = C00011 DIMENSION ANAND» BIWoNDs LEI1O014 O(160) | 
3 000011 aCCouaxiP (aCCucl oo) Bi 
Bi TITS) IPCEVNE.1) @0 TO 200 68,0209 
a 00eels CALL RMIT (oNe}.0Q) 68,0209 
: o0cer) CALL MATOTA® (QeG0N) 68,0209 
c i 
4 200 reo 40303400 pi 
4 0806 40309606 = 
3 wee. meSeasoc a 
i weobonca 40903816 7) 
¥ #1306 60,3808 
90 1 teen mosess2e 
woe GOAC ToL OAC ToL) 0303836 
atliet. : Me909800 
: {lefet n0303850 
BS AO 2 jolett we30386C j 
z POARSP IAT od)? 40303870 } 
wane yo? mesos8ec 
: TF (BoQLLIVAe2ed 40303800 
= 3 ale? 40383600 
: pKtBreJ 0303610 
F 2 CONT rNUE wo3e3e2z0 
7 TF(OcLI ew) Lele 030363¢ 
aa ry wedtr) mesesese 
9 1330 40303650 
00eles i CONTINUE 0303660 
e00107 TPCTTL.L0,0) G0 7O 31 @0,1265 
¢c 
eeell0 HaACCOSERTE (2, 9M) /FLOATF (ND mo30367o 
¢ 
| 3 000123 30) tTeLe (EIT) 40303660 
a 000126 Jselte #0363600 
] 000127 UmAIT Lott) me303700 
000133 yaa(JJeIT) mesosvic 
‘ 000137 ved (JJodJ) mo303720 
rq 000162 wexey 0303736 
if 0003460 Te oSOFOIV( (WeSGQRTF (Wowed, OVOY) )/¥) 
: # 000157 weSQaTr (1, oT#7) e3e3750 ; 
4 ° 000166 aerPOrvitsw) i 
# 000370 reFOtvil.sw) a 
i ® 000176 ccecec 0309780 a 
a qSaSes 40303700 3 
. eCaSece2, 0303800 A 
: 180, 40303016 i 
: 290, 40303620 “a 
e r) W006 mo303030 i 
‘ ‘ eeoze3 wReNRel m03096e0 it 
i 000230 nO 27 fa},N 40303880 j 
4 000233 TMU PoLEIGe Ig 32 mo203e00 
000233 10 WALT TI) wesose70 : 
oee220 vedl yet) ™q309800 yf 
; oee2es ruUegevec meseseee b 
§ eecr27 atl tiee 40203000 : 
a 00233 1F CAMSF (EC) oQLS 1S o1S036 9309010 
4 e002%e 1e = g2eansr (es mesesoze 
eoezee fier mosesese 
4 eoezee 15 PaVOgeurc mosesese 


4-85 


ae ioe 


OT RSLS PR atai_-vtaaheggeaternna trae 


eeerzee alJJo ti OF #939090S0 
00253 TF (ABSE (F) 02) 900016 40363960 ~ 
! e le zedesr (Fr) 9303076 ( 
Es t2e! mo30308e 
- eocze2 ao To @ 40303000 t 
a * ¢ 1 
i. 00262 M2 ACBL Tp wSSexesCeveCCoz mO309000 
: eee27e attedl noseeeie ( i 
ae ocer76 LE CD aE) wesoseze 
: : eeesee a0 To 8 mesceese | 
c ‘ 
i oeese: 12 HCL JJIIP 1G e19 m0300960 (2 
he Vy eA IT) ™0900080 ; 
z vedtJdel) "9904060 aj 
3 deseuscey nosecore 
aT actortiee 40306080 
: ; 1F (ABSF (E) oO11)91501Se22 moseseeo 
mu 21 wR 40364100 
' . (2) eaBSF ce) meseesie 
& a0 To 18 003062320 q 
a 16 ALJ e TICCCONeSCOyeSSoz? 40300130 
b atlotreo. mO3041460 
E 00356 atlbeak 40306150 
& LRT ef2 403046160 
e ao Tn @ M0306) 70 
19 wht yett 40304180 
veAlteJJ) 0304100 
eauegovec woseszeo § 
PavesgaOC ™e304216 
atlortia& me306220 fy 
a(Tesu OF 40306236 Ft 
GeMAx1F (AQSF(E) sABSP (FD) 40306200 : 
TF (600(3)9909013 
13. 9 tlee "0304260 
TF CAASE (£) eABSF (F)) 23024026 03046270 
2e «UK tCS elt 46306280 
z ao to 9 40306200 
000035 23 UK (2) muJ 40304300 
000400 9 TF (Q(T) 0W} 40025025 40306310 
000006 25 waQtry 0304320 
000067 titer 403046330 7 ; 
000050 40 «= TF (LV) 27627033 40304300 ‘ 
000482 330 eB teIT) 40304350 
000457 veB8(f JJ) 40306360 
000063 AlTolL) suegevec 40306370 ia 
00067) a(losJ) eveg~uec 40306380 ey 
000077 27 = CONTINUE 40306300 4 | 
j o0eso2 TP (Wold) 31693030 40300600 i 
7 000506 31 th (1¥)55_S5056 me300410 
ee eoesce $6 AO So TeleN 40304420 
i 000510 100} £20 40304430 
: eoesi2 00 $1 Jaren | MO304460 
: 000513 TP (A( Jo J) 0151051052 ahead 
000820 52 wAlded) ry 
z oeog24 uaJ 49306670 vs 
f oeesas Be) CONTINUE 40300480 
we ° 
ie 000530 TF (Kal) 59050053 40306000 
b e 000832 53 a (Komp@a(fed) we304500 
Z 000se) at3et»eu #0304810 
’ e0eses NO Se Jalen mosessze 
4 000See 19086 eK) #0306830 
4 eeess2 Ql JonreO( Jel) me3008e0 
: e eeesee So = MC do TD SU me30a5se 
e0eses Se CONTINUE uesesses 
¢ ; 
. € ss ruzS CHANGE FORCES FIRST NONAZERO ELEWENT OF EACH EIGENVECTOR PesrTivE, 
89 17 (2¥cL€.0) RETURN 
r 60 00 76 Jw len : 
‘ 0 628 5 @ 10oN { 
JP (O¢8oJ)) Garb2e7e 
62 CONTINUE } 
64 10 64 ft «© 1eN ' 
, ) 66 AlT94) © oBlfeu) 
€ eeeei7 ‘ 70 CONTINUE 


eeeez2 Return 
000ez2 eno =yacost 


VER 10) 


000000 
000001 


000001 
000002 


000003 
000004 


000005 


900006 


000007 
000030 


00091) 
000012 


MacLean 
000012 
00n000 000012 
1301141409220000000)4 


6170090001 
67307 
6140000034 
6150000077 
Sollo 
6325 
7130000000 


63400 
9754000010 « 
67227 
es2orgnucs * 
10666 
aeed 
T5736 
56710 
M20onooori « 


034813 


TOENT ALLER 
PACGRAM LENGTN 


@LOCKS 
PROGRAM® LOCAL 
ENTRY POINTS 


00000) KALLER 


08730/68 


PAGE NO, 


e SOURCE OATE 6820209 CONVERT 70 COC 6e00 
EnTRY RALLER 
Wo JO/OLMALLER D267) 
e 
° WHEN CALLED WITH ARGUMENT Ny RETURNS NAME OF NeTH ROUTINE 
bd Im THE STRING OF CALLING ROUTINESs SUCH THal ZERO 
e CORRESPONDS TO THE ROUTINE THAT CALLED KALLERe 
° ALSO0 RETURNS REL ADOR OF CALL IN AROUMENT 
«e 
KALLER essz + 
$67 1 
$63 87 
S8e 30 
ses 776 
Sai a 
$82 KY 
Sk3 RALLER~] 
Loop $a2 13°87 FETCH yUmPegaCa INSTRUCTION 
Lae ROeKk2 CIRC SHIFT seo RET AOOR NOW IN LOWER ENN 
Sa3 02083 FETCH RETURN@ JUMP INSTRUCTION 
She nQ FETCH Name sNaRGS 
S84 Re NUMBER OF ARGUMENTS 
ut 5 0B4 eDONE 
$82 2-87 OECREMENY COUNTER 
PL A2-Lo0P 
DONE Bks xe PUT NAME TH XO 
SBe ae 
SA? £3-86 FIND @EL a0nR 
Sa? et STOWE IN ARGUMENT N 
we ALLER 
ENN 


JNUSED STORAGE 


22 STATEMENTS 


3. SvMMOLS 


4-87 


we 


SaNENMA on oki 


"ee AST DI RTA = 


Sau ee opt SS 


aye agent 


000006 
eo00ue 
0000U4 
000006 
000008 
000004 
000006 


000008 


ooonen 
000023 
900025 
oo0nas 
000032 
000038 


000102 


000112 
000114 


¢ 
¢c 
¢ 
¢c 
¢ 
¢ 
¢ 
¢ 
c 
¢ 
Cc 
Cc 
¢ 
¢ 
¢C 
¢ 
¢c 
¢c 
c 
¢ 
¢ 
¢ 
¢C 
Cc 
¢ 
€ 
¢ 
Cc 
¢C 
Cc 
Cc 
¢ 
Cc 
¢ 


FUNCTION RORVERG (EMWONe VARRL WORK) 

SuouRce UATE 6960723 SPECIAL HANDLING FON MUGE VALUES OF ENROR 
SOURCE UATE 94920709 ASSURE PRCSN POSITIVEs ADL ENTRY SETVERG 
SOURCE LATE 4660Kia INTULALSZE TO ZEROS 

SOumce VATE 8721208  aDAPTED FROM LIVERGE OF 07.0714 


VAKBL 23 WESETs TO DIVE EMHON TO ZERUcee RONVENG RETURNS A FLAG 
TO INDICATE PROGHESS OF CONVERGENCE> 


SELECTION OF STEP IN WARAL 35 AS FOLLORS ooo 
TF ENOW any ELAST BOTH MUGes MEPEAT INITIAL STEP 
18 GITMER ENOW UM ELAST MUGE (NU) BOTH) s maLVE THE INTERVAL 


36 ROTH EXCH AND ELAST ANE REASONAULEs USe SLOPE THRU LaST 2 PUINTS 


RONVERG PUST we INITIALEZEO AS FOLLUES seoee 
CALL KONVSET (OV ePROSNewOKK) 
ON CALL SETVERG (DVePHCSNeWURK) WREHE cove 
OV 1S INITTIaL INCHKEMENT TO BE USED ON VARGL 
PucSw 3% ALLOWABLE TOLERANCE ON VaReL 
WORK 18 BewWORD WOKK@SPACE USED AY KRONVERG 


09.0709 
60.0709 


WKONVERG RETUWNS FUNCTION VALUE @ 0 IF CONVERGENCE PROCEEDING NORMALLY 


#1 I6 CONVERGENCE COMPLETED 
1 IF FAILURE TO CONVERGE 


TYPICAL USE LS AS FOLLOWS coses 
UIWENSION wORK (8) 
CALL ROAVSET (1000 0KERR» WORK) 
Val00. 
TESTSFCN(V) © GUAL 
IFCKONVENG (TESTeVewORK)) 20010950 
(EmRON MESSAGE) 
(PROCEED) 


EQUIVALENCE (AEUolEG) 

OLMENSTON WORK IO) owAn (8) 

EQUIVALENCE (wHK(1)0 ROUNT Do lWAK(2)9 PRCSN ) 
EQUIVALENCE (WRA(3)9 VNOW Do (WRK(4) 9 VLAST ) 
EQUIVALENCE (WHA(S)9 ENOW Jo(wAK(0)¢ ELAST ) 
EQUIVALENCE (wRK(7)s OV YelwRK(H)e OVLAST ) 
Data (LIMTTEQy) o (KOUNTee}) 


TE CKOUNTeLT00 cORe WONK (DI eh Te0o) CALL CEASE 
® (e@he RONVENG MAS FATLED® OR WAS NOT INITIALIZED ©) 
CALL AMIT(@ewOkR WRK) 
KOUNTSKOUNT 6} 
KFUNCE) 
IF (KOUNT eG oL IME Tea) KEFUNCER 
TF (ROUNT eGTeLIMIT) KFUNCOS 
IF (ROUNT CGT e2 cANDe ABSIFULV(ERROK/ENOW)) oGTole 
e eANDe ERRORPENOW GT .0s 
e eANDe ANS(FOIVIENOW/ELAST)) oGTele) KFUNCOO 
GO TA (20910510010) exFUNC 


10 CALL COUN OUT(A) 
WNITE(G0T1) ROUNTeVLASTSELAST oVNOW cENOW se VANAL CERHOR 


09.0723 


69.709 
69.0709 


69.0723 


Se MRALE AMER TSR aa gay 2 eadiaay Se hata sities nites a 


: 
: 3 1 it 
#4 
: 
; 
; { 5] 
oa 
too 
| 
& oOuUL a? Li FumMal(asne taunts CONVEROING geee vAMLABLE tkheUR / 
r x e cftealion ely ‘ 
x Vib jeekveteo37) D 
OH! wz? Tr TREC ACeL te dt GO TO de ! 
uh y.3 Tare Cl Clio? 
Ore yes Coie eR FRR (] oeORCRAIN UF CALLING PKUGRAMS -.) 
H COP ya? ae OR TR teh ee eyO eS) eRFUNC ( 
2 i 
e annynl cl VUES TS VAuR ! 
Bes a PAALRA burstTstnce 
me Dutra Vim EV ARRL 
: OO0y*eSs bihOwae he 
¥ oMolee Th Cer iwetivetee} 8C TU 99 O9eNTeS | 
3 O00}. IFinmnTele od) GL TO # 69.0783 | 
Oonen TE aes lh wd abel eE LS eANUe AyS(ELadsl) eGTelebh is? GO 10 4] 6920723 eat 
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EVECTORS MATRIX. Tre FIRST COLUMN IS THE VECTOR FOR Tre FIRST 
EIGEN VALUE ETC. SEE JACOBI FOR CALCULATION DETAILS. 


OLMENSTON AUNGN) SEVALUESIND cEV CTORS (NON) 60.0209 


SAVE DIAGONAL ELEMENTS IN EVALUES 
NO } Then 

EVALUES (CP) sellers 

CONT [NUE 


CALCULATE EIGEN vVaLUES 
CALL JACURT (Ast ¥ CTORSeNel Eo oedell) 68.0209 


SwITCH OTAGONAL TERMS AND EIGEN VALUES 
DO 2 Jelen 

WORK ME VALUES (1) 

EVALUES (I) BACT OT) 

ACT,y) se0RK 

CONT TNUE 


MESTORE OFF CIAGONAL TERMS TO a UPPER RIGHT TRIANGLE IS UNCHANGED 
AY JacOBl 

mane} 

90 4 Telem 

nape} 

00 3 selet 

A(Ke cI BA JOR) 

CONT J NUE 

RETURN 

FNC MATIGEN 


MATIGEN 


pod 


cooole 


oo001e 
000014 
000016 
000022 


000023 
o0ones 
000026 
900031 
000033 
000036 
000046 
000050 
000051 
000055 
900060 
090043 


000066 
a0g079 
00007? 
000073 
000077 
000196 
000;10 
00012} 
900°34 
OOO 117 
ONL Ae 
GOO1 
0091 3+ 
00014) 
GO9146 


00038) 
000153 
ong1ss 
Oory re 
oon1?s 
009) 77 


9002716 
000220 
000222 
000233 
0002736 
000236 
g0028) 


oaNe 


one 


nm 


aan) 


Wan 


aon 


oan 


ca) 


60 
54 
aS 


150 
us 


140 


380 


aro 


460 


Suan uTINE MATINV(AeNe NMOS SBomMePTVOTe (PIVOT s TINDER DE TERM) 
SOURCE DATE 620209 CONVERT TO CDC oad0 

£) COUR MATINY MaTHITR INVERSTON 

MaTHTX INWERSTON WITH ACCOMPANYING SOLUTIUN OF LINEAR EdUaTIONS 


MATMTK [MVE ASTON mI Th ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 
SEMEN DTON TRE VUT ON) A UNMAR NMA) oBINMAK G1) gS INDEX (N92) pPIVOTIN 
uf TPard dato 


ETE A Mele 
ye er yey aN 
IPIvaTiodan 
JO SSU Leper 


SEARCH FOR FEVOT ELEMENT 


Awahe0.0 

AD YAS Jey eN 

IE (JOPVOTCUi=b) BOs L0S* 60 
D2 170 Kalen 

TE (LRT ¥OTIRI #1) BOe 1004 7en 
TE CAwSF CAMARA) ARSE (a1 JeK) IE) ASe 100+ 100 
TQUway 

TCOLumen 

AVAKSAL eR) 

CONTINUE 

coNTINUe 

TILVOT EPP ULM) MPO LTVOT YT COLUM) 0) 


"=~ TEWCMANGE WOWS TU PUT PIVOT ELEMENT ON OT AGONAL 


Te {trum TCOLUM) 1402 260¢ 144 
or Tease Ge TEU 
~ enh wma 
se tPaeAd IROweLt 
Se feemen) wal ECOL meLa 
Me shear 
boy cole PLU 
2s Liye 
SwaPani [howe 
she mellem (FCO MeL) 


ae {Ou 
“et 


mis 07 lM yi ) Swab 

re ba Depp elec 

“ane ne deer stocuum 

ervOte pd eatierumes ll COLUM) 


“p TE wma ETE MwOrTVGT CT) 
feloe @lyel wie MY BIVOT FLEMENT 


BTCA me fFULVM) al en 
YG LAs 
AEC Leer PRPL EV COUT CULUMELI ZR IVOT EID) 
Tha gH, Jae WHO 
gre Letom 
ath C CMe PERO EVE CHLUMeLI SHI VOT CED 


HENUCE NONAPTVOT ROWS 


Of 556 Lielen 

TFCLIS1COLUM) 4005 $8Me 000 
Tealedy tCoLum) 

Othe FCOLUM) 0000 

9 460 LayeN 
ACLLeLIMATLU GLI SALICOLUMSL) OT 
TF 1m) $50, SSMe 460 


60000000 


00000010 
40000020 
00000040 
20000020 
00000040 
00000050 
00000060 
00000070 
00000080 


20000110 
00000130 


000001460 
60000189 


00000190 
19000210 
9000220 
0NO%2S0 
90000260 
*00002KR0 
00000290 
0000300 


90000320 


20000 3K0 


00000380 


nao0reic 
en0gcaen 


CON0IAKO 
9000-470 
eoca0aan 


1000NS aA 


0000550 
99000560 
40000570 
40000580 


90000600 


90006630 


MAT INV 


4-107 


o002%4 00 500 Laye™ 06000680 
0002%6 ACL eC SAIL LoL BCI COLUM.) OF 00000660 
000273 CONTINUE 006000670 
00000éuc 

INTERCHANGE COLUMNS 00600690 

00000700 


ES ie een ere 


ano 


000300 9 710 leper 

900301 Lanelel 

000703 TF CINDER CLS TDSTNDER(L02)) 6300 THO 630 

000316 5 jeOwelnvewtl sl) 10000740 
000322 SCILOMBEINDEA LL 02) 

000326 rg 79S nepen 

000327 SwaPeai ks JROM) 

000336 ACK e JRO) BAUR SCOLUM) 

00030) Ame COLUM) &SwaP 

000383 CONT TNue 90000800 


000304 CONTINUE 90000810 
000367 RETURN 29000820 


000350 END maATING 


MATMULT f 


° 
, : 
SUHMUUTINE MATMULT (hes CeNR&sNCARB NCB) i, 
€ YIURCE OATE 6720726 seAND NEW CODE 
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¢ MEMRINGSs Co te GENERAL MEAL MATRIX MULTIPLICATION 
c c =k & Bp WHERE 2p82AaND C ARE MATRICES 
e c NR @ NUMBER OF ROWS IN & CAND C) 
e C NCAMH © NUMGER OF COLUMNS IN & AND ROWS In 9 
¢ NCB ® NUMBER OF COLUMNS IN B (ANN C) 
C 
t ooe NOTE ©0® mMATRIR C MUST WOT OCCUPY THE Same STORAGE AS 
Cc maTHICeS & OR 8 IN THE CALLING SEQUENCE. 
c 
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€ maT P BND P PERMUTES ANO/OM PARTITIONS MATRICES 
@ ¢ 
Py Cc @o15 AN MwA & NCA MATRIX 
Cc 4 1S #& SUR MATRIX UF Ae CHOSEN ACCORDING TO KODER AND KONEC FOR 
Cc WORS AND COLUMNS RESPECTIVELY. 
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c 4OwS AND COLUMNS OF @ MUST CORRESPOND TO THE NUMRER OF DIGITS IN 
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SUBROUTINE PROJ (408.0) 
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¢ 
; c RETURIS VECTOR C = THE PROJECTION OF a ONTO B NIRECTION 
€ 
e 000005 DIMENSTON 0(¢390B 63) 0003) 
¢ 
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000022 
900026 
060030 
0009004 


000068 
000047 
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BRAND NEw ODE 


aRQUMENT ake 3S JONOREN 
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CALL MCNAQIK MESH LOK AR,1) 
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CONTINUE 
NWO8(KOF1 716 
PRINT LetmES@(m) oie) enw) 
FORMAT (Stwo@e 108101 


90 24 
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Levey 


1 
€v=) 
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ic SOURCE OaTE A%e1129 CALL FRASFY WITH NOEAIVY « I 
c S.URCE DaTF 67,0708 REPAIR TRNGEMOING FOR KONROS2 
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ra vECTOR (NOY RELATIVE To RADAR) { 
c 
c eee NUMBER NEGATIVE { 
c oe TRACK 15 RETURNED IN 840aR POLAR COOADINATES, BELATIVE 
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000053 FALL IMTT (begs F ARO) 67,0626 

990085 7 1s N@yeto 20ne 
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SUAROOTINE RORAZ (MONE VALUES S1,S220RBEL} 
¢ SOURCE Date 69.1205 USE KONVERG@ INGTEAD OF AIVERGE 
e SOURCE OATF 62.0209 CONVERT 70 COC 4406 
Cc SOURCE Date 64,1215 SFY UP cai), ALSO 
Cc 
¢ GETURNS ORATTAL ELEMENTS BETWEEN $1 AND $2¢ JUST LIke SUAR ORB? 
ie FxCEor Twar $2 1S TAawEN TO BE ON & ROTATING PLANET. 
Sa fB2OSTTIONS §1 aNC $2 aRE aSSUMEO TO BE CONTEMBORARY) 
vy - POR “ORE wo 1, val UE TS FULGNTST IME 
e v COR MODE @ 2, VAL CE TS EXCESS FLIGHTSTIME (ABOVE MINSENERGY QBIT) 
c FOR VIDE w 3, VALUE 1$ SPEED AT £1 weececesenes INERTLAL 
c FCR “ODE » @, VALUE [$$ VEL*ELEV aNGLE aT $1 4.6 INERTIAL 
ig FOR VIDE w Se VA JE £§ SPEED AT $2 seccvcccceee REL YO BROUND 
¢ FOR VOCE w 6 VALUE 78 VEL@ELEY ANGLE AT $1 ,,. REL TO GROUND 
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) 900359 “OMMON SOANCONS . 
i or, gan, Sive sur, Sue. RBODYs 
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‘ 900010 STMENSTON $3119 492110) eORMEL (10) 
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600a;¢ MEMENSTON 1100493) ONQOVEL (3) oV 13) S40) WORK (BS 66,3298 
eo 
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ae 000083 PALL SITTER (SZ, TIMTRY, $1 me} 
O0Coss VAL @y ALUPF 
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o00C7] «SUNT AS 
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